A6-A120  659  APPLICATION  OF  A  NUMERICAL  SEA  ICE  MODEL  TO  THE  EAST  1/1 

GREENLAND  AREA(U>  COLD  REGIONS  RESEARCH  AND  ENGINEERING 
LAB  HANOVER  NH  U  8  TUCKER  AUG  82  CRRE— 82-16 

F/G  8/12  NL 


UNCLASSIFIED 


1.0  B“ 


m  ...  wOm 

U  In 

B  La  iZO 


12  ui 


11.8 

warn 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS-1963-A 


1.0  Bf“  «** 

=====  «  La  12.2 


E  H4  ■■ 

Sf  k&  izo 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS-1963-A 


in  SJU*  *** 

“  la  |L2 

.  .  «  Itf  |2£ 


L25  ■  1.4  ■  1.6 


1.0  B“  ■*» 

S  w  |2|2 

""  '  ’  E  |j£  ■■ 

.  .  !  N  IM 


jy^  l 


MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  STAN0ARDS-19W-A  r- 


liS'gE 

n  l.l  L“  ■** 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANOAROS-I963-A 


MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  STAN0NWK-19W-A 


■WW/NWANM  mu**.  NLWWM 


WwW 


i  i 


REPORT  82-1 6 

I 

M)A  120  659 


X;.  OCT  2  5  1982 


US  Army  Corps 
of  Engineers 

Cold  Regions  Research  & 
Engineering  Laboratory 


Application  of  a  numerical  sea  ice  model 
to  the  East  Greenland  area 


For  conversion  of  SI  metric  units  to  U.SJ 
British  customary  units  of  measurement 
consult  ASTM  Standard  £380,  Metric  Prac¬ 
tice  Guide,  published  by  the  American  Socie¬ 
ty  for  Testing  and  Materials,  1916  Race  St., 
Philadelphia,  Pa.  19103. 


Cover:  Vector  representation  of  average 
wind  field. 


CRREL  Report  82-16 

August  1982 


Application  of  a  numerical  sea  ice  model 
to  the  East  Greenland  area 

Walter  B.  Tucker,  III 


Prepared  for 

OFFICE  OF  NAVAL  RESEARCH 
Approved  for  public  release;  distribution  unlimited. 


i  innasatned 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


(  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS  jl 

BEFORE  COMPLETING  FORM  I 

4.  TITLE  fend  Subtitle)  r 

APPLICATION  OF  A  NUMERICAL  SEA  ICE  MODEL 

TO  THE  EAST  GREENLAND  AREA 

S.  TYPE  OF  REPORT  4  PERIOD  COVERED 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORS 

Walter  B.  Tucker,  III 

S.  CONTRACT  OR  GRANT  NUMBER!*) 

t.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 

Hanover,  New  Hampshire  03755 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 

It-  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 

Arlington,  Virginia  22217 

12.  REPORT  DATE 

August  1982 

IS.  NUMBER  OF  PAGES 

49 

14.  MONITORING  AGENCY  NAME  4  AODRESSfff  dltloront  from  Controlling  OU!  cm) 

IS.  SECURITY  CLASS,  (ol  thte  report) 

Unclassified 

j^j|||j.wui4A;M:ir,Ti  m  i  ii  i  m l 

l«.  DISTRIBUTION  STATEMENT  (ol  thin  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  (ha  abafract  an  tar*  d  In  Black  39,  It  dWefmrt  from  Report) 

IS.  SUPPLEMENTARY  NOTES 

It-  KEY  WORDS  (Continue  on  reeeree  a/da  if  neceeeetr  and  Identify  by  block  number) 

Greenland 

Ice 

Models 

Ocean  models 

Sea  ice 

30.  ABSTRACT  fOelbia  am  rarer aa  afrk  H  neeeeearr  and  Identify  bjr  block  number) 

*  A  dynamic-thermodynamic  sea  ice  model  which  employs  a  viscous-plastic  constitutive  law  has  been  applied  to  the 

East  Greenland  area.  The  model  is  run  on  a  40-km  spatial  scale  at  1  /4-day  time  steps  for  a  60-day  period  with  forcing 
data  beginning  on  1  October  1979.  Results  tend  to  verify  that  the  model  predicts  reasonable  thicluiesses  and  velocities 
within  the  ice  margin.  Thermodynamic  ice  growth  produces  excessive  ice  extent,  however,  probably  due  to  inadequate 
parameterization  of  oceanic  heat  flux.  Ice  velocities  near  the  free  ice  edge  are  also  not  well  simulated,  and  preliminary 
investigations  attribute  this  to  an  improper  wind  field  in  this  area.  A  simulation  which  neglects  ice  strength,  effectively 
damping  ice  interaction  with  itself  and  allowing  no  resistance  to  deformation,  produces  excessive  ice  drift  toward  the 
coast  and  results  in  unrealistic  nearshore  thicknesses.  A  dynamics-only  simulation  produced  reasonable  results,  including  - 

DO  «  jnTti  1473  EornoN  of  *  now  is  is  obsolete 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWfcen  Data  Bn farad) 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfRhan  Data  Bntand) 


Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGEflVhwt  Data  Bntarat t) 


PREFACE 

This  report  was  prepared  by  Walter  B.  Tucker,  111,  Geologist,  Snow  and  Ice 
Branch,  Research  Division,  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory.  The  work  was  supported  by  the  Arctic  Program  of  the  Office  of  Naval 
Research. 

The  author  is  indebted  to  Dr.  W.D.  Hibler,  III,  for  providing  the  computer  code, 
ice  growth  rates,  and  much  guidance  in  the  early  stages  of  this  work.  The  guidance, 
encouragement  and  comments  of  Professors  R.L.  Haney  and  R.G.  Paquette  of  the 
Naval  Postgraduate  School,  Monterey,  California,  are  also  gratefully  acknowledged. 
The  author  also  thanks  E.  Perkins  and  M.  Pacillo  for  much  time-consuming  illustra¬ 
tion  preparation,  and  especially  D.  Harp  for  an  excellent  job  of  typing  the  original 
manuscript. 


CONTENTS 


Abstract . 

Preface . 

Introduction . 

Model  description  and  application . 

Results  and  discussion . 

Wind  and  current  fields . 

Standard  simulation . 

Thermodynamic  simulation . 

Zero  ice  strength . 

Zero  ice  import . 

Zero  currents . 

Modified  currents . 

Zero  winds . : . . 

Dynamics  simulation . 

Summary  and  concluding  remarks. 
Literature  cited . 


ILLUSTRATIONS 


Figure 

1 .  Model  grid  with  solid  and  free  boundaries . 

2.  Model  flow . 

3.  Wind  and  current  fields . 

4.  Late  summer  currents,  average  ice  margins,  and  manned  ice  station 
ck 


5.  Simulated  ice  thickness  fields . 

6.  Simulated  ice  compactness  fields . 

7.  Averaged  simulated  ice  velocities . 

8.  Five-day  averaged  velocities  for  the  row  of  grid  points  immediately 

above  the  1 1 .0-m  coastal  ice  buildup  shown  in  Figure  Sg . 

9.  Five-day  averaged  velocities  for  the  row  of  grid  points  immediately 

below  the  1 1 .0-m  coastal  ice  buildup  shown  in  Figure  5g . 

10.  Simulated  daily  changes  in  total  ice  volume  and  daily  volumes  of 

growth,  northern  inflow  and  southern  outflow . 

11.  Trajectories  of  ICEX  buoys  1564  and  1568 . 

12.  Observed  and  predicted  components  of  velocity  for  ICEX  buoy  1564. . 

13.  Observed  and  predicted  components  of  velocity  for  ICEX  buoy  1568. . 

14.  Actual  and  predicted  trajectories  for  buoys  1564  and  1568 . 

15.  Actual  and  predicted  trajectories  for  buoy  1568 . 

16.  Actual  and  cumulative  daily  predicted  trajectories  for  buoys  1564  and 

1568  . . . 

17.  Actual  buoy  trajectories  superimposed  on  the  60-day  averaged  ice  vel¬ 
ocity  field . 


m 

u 


i 


Figure 

18.  Sixty-day  thickness  and  compactness  fields  for  the  thermodynamic  sim¬ 
ulation  . 

19.  Sixty-day  thickness  and  compactness  fields  for  the  zero-strength  simu¬ 
lation  . . .  .• . 

20.  Sixty-day  averaged  velocity  field  for  the  zero-strength  simulation . 

21 .  Actual  and  cumulative  daily  predicted  trajectories  for  buoy  1564  for  the 

zero-strength  simulation .  . 

22.  Sixty-day  thickness  and  compactness  fields  for  the  zero-inflow  simula¬ 
tion  . 

23.  Sixty-day  averaged  velocity  field  for  the  zero-inflow  simulation . 

24.  Actual  and  cumulative  daily  predicted  trajectories  for  buoy  1564  for  the 

zero-inflow  simulation . 

25.  Sixty-day  averaged  ice  velocities  for  the  zero-current  simulation . 

26.  Actual  and  cumulative  daily  predicted  trajectories  for  buoy  1564  for  the 

zero-currents  simulation . 

27.  Sixty-day  averaged  velocity  field  for  the  modified-currents  simulation. 

28.  Sixty-day  thickness  and  compactness  fields  for  the  modified-currents 

simulation . 

29.  Actual  and  cumulative  daily  predicted  trajectories  for  buoy  1564  for  the 

modified  currents  simulation . 

30.  Sixty-day  thickness  and  compactness  fields  for  the  zero-winds  simula¬ 
tion  . 

31.  Sixty-day  averaged  velocity  field  for  the  zero-winds  simulation . 

32.  Actual  and  cumulative  daily  predicted  trajectories  for  buoy  1564  for  the 

zero-winds  simulation . 

33.  Sixty-day  thickness  and  compactness  fields  for  the  dynamics  simulation 

34.  Sixty-day  averaged  velocity  field  for  the  dynamics  simulation . 

35.  Actual  and  cumulative  daily  predicted  trajectories  for  buoy  1564  for  the 

dynamics  simulation . 


TABLES 

Table 

1.  Areas  of  ice  cover  for  the  standard  simulation . 

2.  Areas  of  ice  cover  for  the  thermodynamic  simulation . 

3.  Areas  of  ice  cover  for  the  dynamics  simulation . 

4.  Model  sensitivity  test  mass  balance  results . 

5.  Simulated  velocity  comparisons  with  buoy  1564 . 


APPLICATION  OF  A  NUMERICAL  SEA  ICE  MODEL 
TO  THE  EAST  GREENLAND  AREA 


Walter  B.  Tucker,  III 


INTRODUCTION 

The  Greenland  Sea  is  an  area  of  confluence  for 
polar  and  temperate  systems  of  both  the  atmos¬ 
phere  and  ocean.  In  the  atmosphere,  migratory 
cyclones  formed  over  or  near  North  America  track 
through  the  Greenland  and  Norwegian  Seas,  fre¬ 
quently  undergoing  cyclolysis.  These  dying 
cyclones  give  rise  to  the  Icelandic  Low  (Sanders 
and  Gyakum  1980).  In  the  hydrosphere,  a  com¬ 
plex  system  of  currents  results  when  the  warm 
North  Atlantic  current,  a  Gulf  Stream  outflow 
that  flows  into  the  Arctic  Basin  west  of  Spits¬ 
bergen,  meets  the  cold  East  Greenland  current 
that  flows  south  out  of  the  Arctic  Basin  along  the 
Greenland  coast. 

Along  with  the  colder,  less  dense  water  that  is 
transported  into  the  region  by  the  East  Greenland 
current,  sea  ice  is  advected  out  of  the  Arctic  Basin. 
This  ice  transport  is  greatly  assisted  and  possibly 
dominated  by  the  generally  northerly  winds  which 
result  from  the  cyclones  transiting  and  stagnating 
over  the  Greenland  and  Norwegian  Seas.  In  addi¬ 
tion,  the  ocean  surface  heat  balance  is  favorable 
for  the  production  of  new  ice  here  during  the  win¬ 
ter  months,  a  factor  which  further  increases  the  ice 
extent  during  this  period.  Summer  warming  ab¬ 
lates  the  ice  being  transported  south,  resulting  in  a 
greatly  reduced  ice  extent.  The  presence  of  this 
seasonally  varying  sea  ice  cover  results  in  a  system 
of  highly  complex  air-sea  interactions  and  feed¬ 
back  effects  that  are  not  well  understood. 

The  East  Greenland  area  is  of  interest  to  many 
nations  for  economic  and  military  reasons.  Sea  ice 


severely  hampers  surface  navigation,  affecting 
both  the  commercial  maritime  and  fishing  indus¬ 
tries.  Subsurface  navigation  is  affected  by  the  in¬ 
fluence  of  sea  ice  on  the  acoustic  regime.  Hydro¬ 
acoustic  instrumentation  is  affected  by  ice  scatter¬ 
ing  and  high  ambient  noise  levels  due  to  ice  floe 
collisions  and  deformation  (Kozo  and  Diachok 
1973,  Diachok  and  Winokur  1974). 

It  is  obvious  that  there  is  a  need  to  understand 
the  effects  of  sea  ice  on  oceanic  and  atmospheric 
processes.  A  logical  first  step  is  to  attempt  to 
understand  which  processes  control  the  presence 
and  variability  of  the  sea  ice  in  this  region.  With 
this  understanding,  atmospheric  and  oceanic 
models  could  be  improved  by  including  the  re¬ 
sponse  of  the  ice  cover  to  predicted  results,  there¬ 
by  implementing  crude  feedback  mechanisms. 
Later,  coupling  of  these  models  would  begin  to 
delineate  the  more  complex  processes. 

The  major  components  which  govern  the  sea  ice 
balance  in  any  ice-covered  region  are  the  thermo¬ 
dynamic  balance  at  the  sea  surface,  the  air  and  wa¬ 
ter  stresses  upon  the  ice,  the  Coriolis  force,  and 
the  internal  ice  stress  (the  stress  transmitted  by  the 
ice  itself)-  Consideration  of  the  East  Greenland 
area  as  a  separate  entity  also  requires  that  the  flux 
of  ice  into  the  region  from  the  Arctic  Basin  be  in¬ 
cluded  as  a  component  in  the  mass  balance  of  ice. 
The  role  of  each  of  these  components  with  respect 
to  the  East  Greenland  region  has  not  been  made 
clear. 

Previo-  is  studies  of  tf  ice  balance  in  this  region 
have  foe  "rima  on  variations  in  ice  extent. 
Both  season  anu  .nterannual  variations  have 
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been  examined  in  these  predominantly  climatolo¬ 
gical  studies.  Vowinckel  (1964)  concluded  that  the 
international  variations  in  ice  extent  are  less  than 
seasonal  variations.  He  reasoned  that  potentially 
large  year-to-year  variations  are  generally  counter¬ 
balanced  by  extremes  in  freezing  and  thawing. 
That  is,  the  greater  the  extent  in  April,  the  greater 
the  reduction  of  ice  by  melting  in  the  summer.  His 
study  also  concluded  that,  on  the  average,  season¬ 
al  variations  are  due  to  fluctuations  in  the  amount 
of  ice  being  imported  from  the  Arctic  Basin.  He 
estimated  the  total  southward  ice  transport  for 
each  month  by  examining  ice  extent  charts  and  ap¬ 
plying  simple  assumptions  concerning  freezing 
and  thawing.  The  transport  by  wind  alone  was 
then  estimated  by  applying  Zubov’s  formula  (Zu¬ 
bov  194S)  to  monthly  pressure  differences  at  cer¬ 
tain  latitudes.  The  transport  attributable  to  cur¬ 
rents  was  taken  to  be  the  difference  between  the 
total  and  wind  transports.  These  calculations 
showed  that  currents  dominated  ice  transport  dur¬ 
ing  the  winter  (September-April).  Wind  transport 
was  approximately  half  of  current  transport  dur¬ 
ing  that  period. 

Skov  (1970)  believes  that  year-to-year  variations 
in  ice  extent  are  caused  by  ocean  current  varia¬ 
tions.  In  particular,  the  fluctuations  in  the  flow  of 
Polar  and  Atlantic  water  into  the  Greenland  Sea 
cause  the  northward  oceanic  heat  transport  to 
vary,  thus  influencing  the  ice  extent.  Aagaard 
(1972),  on  the  other  hand,  found  that  severe  ice 
years  were  accompanied  by  anomalous  atmos¬ 
pheric  pressure  fields.  Using  Sverdrup  dynamics 
he  showed  that  the  southward  transport  of  polar 
water  increased  during  these  anomalous  years, 
bringing  unusually  large  amounts  of  ice  south. 
With  this  reasoning  he  has  attributed  anomalous 
ice  extents  to  variations  in  the  mean  wind  and  cur¬ 
rent  fields  by  assuming  that  the  ocean  will  respond 
within  a  reasonable  time  (several  months)  to  the 
mean  wind  stress  field.  His  hypothesis  is,  more¬ 
over,  based  on  local  forcing. 

Other  large-scale  climatic  studies  have  related 
the  ice  extent  to  various  atmospheric  parameters. 
Walsh  and  Johnson  (1979)  designed  a  study  to  as¬ 
sess  interactions  of  the  sea  ice  and  the  atmosphere 
by  cross  correlating  meteorological  and  ice  extent 
fields  which  were  represented  by  empirical  ortho¬ 
gonal  functions.  In  examining  these  first-order 
feedback  effects,  they  found  that  the  ice  extent  re¬ 
sponded  more  strongly  to  atmospheric  forcing  oc¬ 
curring  one  to  two  months  previously  than  did  the 
atmosphere  to  the  ice  extent  for  any  lag  or  lead 
time  during  an  ice  retreat  period  (summer).  Dur¬ 
ing  the  ice  advance  period,  however,  the  best  cor¬ 


relations  occurred  with  zero  lag,  indicating  iir  me¬ 
diate  forcing,  and  atmospheric  response  to  ice  ex¬ 
tent  was  equally  as  strong  as  the  response  of  the 
ice  extent  to  the  atmosphere.  The  meteorological 
variable  showing  the  highest  correlation  with  ice 
extent  was  generally  the  surface  temperature  field. 
In  a  different  study,  which  examined  specific  me¬ 
teorological  features,  Kelly  (1978)  suggested  that 
ice  extent  in  the  Greenland  Sea  may  be  closely  re¬ 
lated  to  the  position  of  the  Icelandic  Low. 

That  both  local  winds  and  ice  production  play  a 
major  role  in  ice  extent  was  promoted  by  Einars- 
son  (1972).  He  pointed  out  that  since  drifting  ice 
stations  often  moved  across  the  assumed  strong 
current  along  the  Belgica  Bank  of  northeast 
Greenland,  the  wind  must  be  significantly  influ¬ 
encing  the  ice  drift.  In  calculating  an  annual  ice 
mass  budget,  he  found  that  export  out  of  the  Den¬ 
mark  Strait  exceeded  the  inflow  at  76  °N.  Thus  the 
region  must  be  a  net  producer  of  ice.  He  also 
agreed  with  Vowinckel  (1964)  in  finding  that  win¬ 
ter  drift  rates  are  much  greater  than  those  of  other 
seasons  due  to  a  stronger  northeasterly  pressure 
gradient. 

Short-term  rapid  advances  of  the  ice  edge  have 
been  attributed  to  different  mechanisms.  For  in¬ 
stance,  Einarsson  (1972)  cited  specific  studies 
which  have  found  rapid  advances  of  ice  in  the  vi¬ 
cinity  of  Iceland  to  be  preceded  by  tongues  of  low 
salinity  water.  A  rapid-advance  feature  has  been 
noted  for  centuries  by  fishermen  in  the  vicinity  of 
Jan  Mayen.  A  large  cape-like  extension  of  ice  (apt¬ 
ly  named  Odden)  protrudes  northeastward,  delin¬ 
eating  a  bay  of  open  water  to  the  northwest  (Nord- 
bukta).  Vinje  (1977)  believes  that  the  causes  of  this 
rapid  advance  are  a  weakened  Icelandic  Low 
(which  would  lessen  the  easterly  winds)  and  a  well- 
developed  oceanic  circulation  (which  would  trans¬ 
port  the  ice  to  the  east).  In  a  separate  study  San¬ 
derson  (1971)  found  that  rapid  advances  of  the  ice 
edge  in  the  Greenland  Sea  were  not  accounted  for 
by  monthly  mean  winds,  currents  or  ice  growth 
rates.  Instead,  he  found  a  significant  correlation 
between  ice  edge  advance  and  the  monthly  mean 
wind  anomaly  from  the  northwest  quadrant.  This 
anomaly,  which  constitutes  the  departure  of  the 
wind  from  its  normal  value,  presumably  boosts 
the  southeasterly  branch  of  the  East  Greenland 
Current,  thus  stimulating  a  large  ice  transport  to 
the  east. 

In  addition  to  these  synoptic  scale  processes 
which  contribute  to  ice  drift  and  extent,  Wadhams 
(1980a)  points  out  that  smaller  scale  processes  may 
contribute  significantly  to  the  ice  extent.  He  be¬ 
lieves  that  wave-induced  pulverization  of  ice  near 


n  conjunction  with  an  off-ice  wind 
e  the  pulverized  ice  to  melt  very  rapid- 
*cess  may  cause  a  large  enough  effect  to 
irameterization  in  a  model  which  pre- 
:  edge  location.  He  also  points  out  that 
iddies  associated  with  an  unstable  oce- 
Front  may  cause  rapid  disintegration  of 
elation  proceeds  by  the  eddies  drawing 
he  warmer  ocean  where  melting  is  sig- 
d  by  moving  floes  into  the  proximity  of 
n  where  pulverization  and  subsequent 
:e  place. 

le  above  investigations  have  found  re¬ 
indicate  that  all  major  components  in 
ince  are  important.  The  question  of  re- 
irtance  of  the  terms  remains,  however. 
;re  are  limits  to  what  can  be  resolved  by 
tudies  in  this  region  (due  primarily  to 
and  low  accuracy  of  observational 
emed  that  an  ice  modeling  study  might 
The  idea  of  this  type  of  investigation 
)  attempt  to  sort  out  the  major  process- 
a  series  of  model  sensitivity  tests. 
(1969)  formulated  the  physical  frame- 
l  sea  ice  model  applicable  to  the  East 
Sea.  This  theoretical  framework  con- 
terms  in  the  momentum  balance  and 
ice  as  an  isotropic  elastic  medium  to 
internal  ice  stresses.  In  addition,  a  con¬ 
ation  accounted  for  ice  concentration 
d  for  growth  and  ablation.  Unfortu- 
model  was  never  taken  beyond  the  for- 
tage.  In  addition  to  Karlsson’s  model- 
several  Russian  investigators  have  ap- 
e  ice  balance  models  to  the  East  Green- 
(Lebedev  and  Uralov  1976,  Antropova 
1977).  These  models  basically  estimat- 
>r  components  affecting  the  ice  balance 
tflow  and  growth)  but  without  a  proper 
)f  the  actual  ice  dynamics  within  the 

>rt  presents  the  results  of  the  applica- 
ynamic-thermodynamic  sea  ice  model 
Greenland  area.  Preliminary  results  of 
ive  previously  been  reported  by  TuckeT 
(1981).  Further  analyses  and  the  results 
tensitivity  tests  are  reported  here.  This 
sents  the  First  attempt  to  apply  a  com- 
!  model  specifically  to  the  East  Green- 


developed  by  Hibler  (1979).  This  particular  model 
was  selected  for  use  in  this  study  because  it  has 
previously  yielded  very  reasonable  results  in  Arctic 
Basin  studies  (Hibler  1979)  and  because  the  nu¬ 
merical  code  has  been  documented  (Hibler  1980a) 
and  can  be  applied  to  any  specific  region  with  rela¬ 
tive  ease. 

Basic  components  of  the  model  include  a  mo¬ 
mentum  balance,  a  constitutive  law,  an  ice  thick¬ 
ness  distribution,  an  ice  strength  parameterization 
and  a  thermodynamic  balance. 

The  momentum  equation  for  ice  floating  on  an 
ocean  is 

m  ^  =  C  +  rw  +  ra  +  F  +  G  +  T  (1) 

where  u  is  the  ice  velocity,  m  is  the  ice  mass  per 
unit  area,  C  is  the  Coriolis  force,  tw  and  Ta  are  the 
water  and  air  stresses,  F  is  the  force  due  to  internal 
ice  stress  variation,  G  is  the  force  due  to  long-term 
geostrophic  currents,  and  T  is  the  force  attributed 
to  the  tilt  of  the  ocean  surface.  The  acceleration 
term  ( m  du/dt),  a  total  derivative,  is  further  bro¬ 
ken  into  the  local  acceleration  plus  momentum  ad- 
vection. 

The  constitutive  law  is  of  the  form 

°ij  (2) 

where  is  a  two-dimensional  stress  tensor,  ejj,  is 
the  strain  rate  tensor,  and  P  is  a  pressure  term  re¬ 
presenting  ice  strength,  which  depends  upon  the 
ice  thickness  distribution.  £  and  17  are  nonlinear 
shear  and  bulk  viscosities,  and  their  values  depend 
on  £j:  and  P  in  accordance  with  a  viscous-plastic 
rheology,  The  details  of  this  constitutive  law  are 
presented  by  Hibler  (1979).  The  law  in  this  form 
allows  the  ice  to  deform  as  a  linear  viscous  (New¬ 
tonian)  fluid  at  small  strain  rates  but  yields  as  a 
purely  plastic  material  at  higher  strain  rates.  The 
usual  or  normal  range  of  strain  rates  causes  fre¬ 
quent  plastic  yielding  as  manifested  by  pressure 
ridge  and  lead  formation.  Once  the  stress  tensor  is 
obtained  from  the  constitutive  relationship,  the 
force  components  due  to  internal  ice  stress  are  cal¬ 
culated  from 

Fj  =  doj/dxj.  (3) 


ESCRIPTION  AND  APPLICATION 

:e  model  utilized  in  this  study  is  a  two- 
vel,  viscous-plastic  model  which  was 


The  ice  strength  parameterization  couples  the 
ice  strength  to  the  thickness  distribution.  The  ice 
strength  pressure  term  P  in  eq  2  is  a  function  of 
thickness  and  compactness  (concentration)  ac¬ 
cording  to 
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P  =  P* /» exp  l -C(l -/!)]. 


(4) 


Here  P*  and  C  are  fixed  empirical  constants,  h  is 
the  average  ice  thickness  for  the  grid  cell,  and  A  is 
the  compactness,  which  represents  the  fractional 
area  of  the  grid  cell  (varying  from  0.0  to  1 .0)  cov¬ 
ered  by  ice  of  thickness  h. 

The  evolution  of  ice  thickness  and  compactness 
is  governed  by  two  continuity  equations: 


dh 
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+  Sh  +  diffusion 


+  SA  + diffusion 


(5) 

(6) 


where  u  and  v  are  velocity  components  in  the  x 
and  y  directions  and  Sh  and  SAare  thermodynamic 
let  ms  which  govern  the  ice  thickness  and  concen¬ 
tration  due  to  growth  and  decay.  The  diffusion 
terms  are  necessary  for  numerical  stability.  The 
thickness  and  compactness  of  ice  in  each  grid  cell 
are  determined  by  eq  5  and  6  for  each  time  step. 
The  remainder  of  the  grid  cell  (fractionally,  1  -  A) 
is  considered  to  be  open  water. 

A  surface  heat  balance  equation,  together  with 
a  simple  thermodynamic  ice  model,  were  used  to 
calculate  the  growth  rates,  Sh,  in  the  manner  de¬ 
scribed  by  Hibler  (1980b).  The  balance  equation 
included  terms  for  incoming  long-  and  shortwave 
radiation,  outgoing  longwave  radiation,  sensible 
and  latent  heat  fluxes,  and  ice  conductivity.  The 
external  data  required  to  solve  this  equation  came 
from  the  National  Climatological  Center’s  (NCC) 
daily  analyzed  fields  (temperature,  humidity  and 
pressure)  and  from  climatological  estimates 
(cloudiness).  Radiation  values  were  calculated  as 
described  by  Hibler  (1980b).  In  a  separate  run,  the 
ice  growth  rates  were  calculated  for  each  grid 
point  at  O.S-m  thickness  levels  (Hibler,  pers. 
comm.),  then  stored  for  later  access  by  the  model 
which  interpolated  a  growth  rate  to  the  proper 
thickness  level.  The  change  in  compactness  due  to 
growth  and  decay,  SA,  is  calculated  as  detailed  by 
Hibler  (1979).  This  effect  is  parameterized  so  as  to 
allow  the  amount  of  open  water  (or  very  thin  ice) 
to  rapidly  decrease  under  growth  conditions  and 
to  slowly  increase  during  periods  of  melting. 

Initial  simulation  runs  which  tested  only  the 
thermodynamic  portion  of  the  numerical  code 
found  ice  growth  to  be  excessive.  This  was  pre¬ 
sumably  due  to  the  lack  of  oceanic  heat  flux.  Wa¬ 
ter  temperature  in  the  heat  balance  equation  was 
specified  to  be  271.2  K,  the  freezing  point  of  sea¬ 
water.  This  implies  that  no  ice  ablation  as  a  result 


of  either  boundary  layer  heat  storage  or  advection 
of  warmer  waters  into  this  area  occurs.  In  light  of 
this,  a  crude  oceanic  heat  flux  was  incorporated  by 
adding  a  0. 1  -m-per-day  decay  rate  to  the  thin  ice 
growth  rates  east  of  a  fixed  boundary  in  the  grid 
area.  West  of  this  boundary  this  melt  rate  supple¬ 
ment  falls  off  in  a  cosine  manner  until,  tinally,  no 
modification  is  made  to  the  growth  rates.  I  he  idea 
of  establishing  this  zone  was  to  attempt  to  crudely 
simulate  the  oceanic  Polar  Front,  which  roughly 
follows  the  ice  margin,  having  relatively  warm  wa¬ 
ters  to  the  east  and  below-freezing  temperatures  to 
the  west  (Wadhams  et  al.  1979,  Aagaard  and 
Coachman  1968).  The  position  of  this  zone  within 
the  model  grid  is  shown  in  Figure  1 .  It  was  posi¬ 
tioned  roughly  according  to  the  location  of  the  ice 
edge  at  the  end  of  November  1979. 

The  air  and  water  stresses  as  shown  in  the  mo¬ 
mentum  equation  (eq  1)  were  calculated  from  sim¬ 
ple  nonlinear  drag  laws  which  assumed  constant 
turning  angles  and  required  geostrophic  winds  and 
ocean  currents.  Winds  were  calculated  from  NCC 
daily  analyzed  sea  level  pressure  fields.  Currents, 
as  well  as  ocean  tilt,  were  determined  from  a  tem¬ 
porally  constant  dynamic  height  field  (Syd  Levi- 
tus,  pers.  comm.).  These  fields  are  discussed  fur¬ 
ther  in  the  Results  and  Discussion  section. 

The  overall  flow  of  the  model  is  shown  in  Figure 
2.  The  simultaneous  equations  1,  4,  5  and  6  are 
solved  by  finite  difference  techniques  using  a  stag¬ 
gered  grid  procedure.  The  momentum  equation  1 
is  solved  by  a  semi-implicit  predictor-corrector 
technique,  with  velocities  being  calculated  by 
over-relaxation  for  each  of  the  two  time  steps.  The 
thickness  and  compactness  continuity  equations 
are  solved  explicitly  with  a  modified  Euler  step. 

A  40-km,  31  x45  grid  covering  the  East  Green¬ 
land  area  was  established  for  the  simulations.  A 
Lambert  azimuthal  equivalent  projection  pro¬ 
vided  grid  cells  of  equal  areas.  The  location  of  the 
grid  with  boundaries  is  shown  in  Figure  1.  Using 
the  staggered  grid  procedure,  velocities  are  calcu¬ 
lated  for  the  grid  points,  with  thicknesses  and 
compactnesses  specified  for  the  cells  between  the 
grid  points. 

The  Fram  Strait  (northwest),  the  Denmark 
Strait  (south)  and  the  entire  eastern  boundary  are 
designated  as  open  boundaries  (through  which  in¬ 
flow  and  outflow  can  take  place).  For  ice  strength 
calculations,  the  thickness  of  the  open  boundary 
cells  in  the  Denmark  Strait  and  the  eastern  boun-' 
dary  is  taken  to  be  the  average  thickness  of  adja¬ 
cent  cells  located  inside  the  boundary.'  For  the 
Fram  Strait,  an  inflow  region,  a  different  proce¬ 
dure  was  followed.  Because  this  area  constitutes 
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Figure  /.  Model  grid  with  solid  and  free  boundaries  (shaded}.  Dashed  lines  represent  bounda¬ 
ries  of  the  oceanic  heat  flux  parameterization. 


the  major  outflow  region  for  the  Arctic  Basin,  the 
thicknesses  of  these  cells  were  specified  indepen¬ 
dent  of  time. 

Velocity  values  at  the  solid  boundaries  are  set  to 
zero.  With  the  viscous-plastic  rheology,  which 
contains  terms  to  account  for  compressive  and 
shear  stress,  effects  of  the  coastline  should  be  ade¬ 
quately  reproduced.  That  is,  longshore  shear  and 
high  strength  areas  due  to  ice  thickness  buildup 
presumably  are  inherent  in  the  plastic  rheology 
without  further  boundary  stipulations.  This  was 
shown  to  be  the  case  in  an  Arctic  Basin  study  with 
this  model  (Hibler  1979). 

The  relatively  small  grid  size  (40  km)  and  large 
ice  velocities  in  this  region  required  the  model  to 
be  run  at  a  '4 -day  (21,600-s)  time  step  to  satisfy 
the  Courant-Fredrichs-Levy  stability  criterion  (in 
this  case  A/<Ax[2(u]  +  vJ)]'l/j.  Winds  and  growth 


rates  were  interpolated  from  their  daily  values  to 
this  interval.  All  other  model  prameters  were  iden¬ 
tical  to  those  used  by  Hibler  (1979)  for  the  Arctic 
Basin  study  with  the  exception  of  the  Coriolis 
parameter  and  P*.  For  this  study,  the  Coriolis 
parameter  was  calculated  for  each  grid  point  loca¬ 
tion  and,  therefore,  it  varies  with  latitude.  The 
constant  P*,  used  for  the  determination  of  ice 
strength  (eq  4),  was  set  to  four  times  the  value  us¬ 
ed  in  the  Arctic  Basin  simulations  (5.0  x  10*  N 
nr').  This  change  was  implemented  when  initial 
tests  showed  ice  velocities  to  be  excessive,  presum¬ 
ably  due  to  the  large  magnitudes  and  variability  of 
the  daily  wind  fields.  The  previous  simulation  had 
used  8-day  averaged  winds,  which  inherently  pro¬ 
vided  spatially  and  temporally  smoothed  fields. 

Because  CRREL’s  computer  was  small,  the  sim¬ 
ulations  were  restricted  to  a  60-day  period  in  order 
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C  Coriolis  Force 
rw  Water  stress  due  to  ice  motion 
rA  Air  stress 

F  Internal  ice  stress  variation 
G  Ocean  currents 
T  Ocean  tilt 
Du 

— 1  Ice  acceleration  and  momentum  advection 
m  ice  mass  per  unit  area 

Figure  2.  Model  flow  (from  Hibler  1980a). 


to  obtain  a  reasonable  turnaround  time.  This  fact 
limited  the  scope  of  this  investigation  to  short¬ 
term  or  seasonal  effects.  The  continuous  period  of 
October  through  November  1979  was  chosen  for 
study,  primarily  because  it  was  a  period  of  rapid 
ice  expansion  and  thus  would  allow  an  assessment 
of  the  relative  importance  of  dynamics  and  ther¬ 


modynamics  to  the  ice  expansion.  In  addition,  po¬ 
sition  data  for  drifting  buoys  located  on  the  ice 
were  available  for  this  time  period  (Kloster  and 
Rafto  1980).  initial  ice  compactness  was  digitized 
from  the  2  October  1979  ict  chart  as  published  by 
the  Naval  Polar  Oceanography  Center  (NPOC 
1979).  Thickness  for  the  initial  field  was  estimated 
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a.  60-day  averaged  geostrophic  wind  field. 


b.  Geostrophic  ocean  currents. 


Figure  3.  Wind  and  current  fields. 


by  allowing  it  to  vary  linearly  with  latitude,  1 .0  m 
at  67  °N  to  3.2  m  at  83  °N.  These  estimates  seemed 
reasonable  based  on  data  reported  from  subma¬ 
rine  transects  of  the  area  during  different  time  pe¬ 
riods  (Kozo  and  Tucker  1974,  Wadhams  1980b). 
Similarly,  thickness  for  the  Fram  Strait  inflow 
cells,  which  remained  constant  for  the  simula¬ 
tions,  was  specified  to  be  3.2  m  for  the  cells  near¬ 
est  the  coast  and  decreased  linearly  to  0.0  m  for 
the  most  northeasterly  cell. 

Various  simulations  were  carried  out  to  assess 
the  response  of  the  model  to  different  forcing  pro¬ 
cesses  during  this  60-day  period.  The  primary  sim¬ 
ulation,  referred  to  as  the  standard  simulation  (or 
run),  incorporates  the  entire  dynamic-thermody¬ 
namic  model  with  forcing  fields  as  described.  Oth¬ 
er  simulations  test  the  sensitivity  of  the  model  1 )  to 
thermodynamics  alone,  2)  to  zero  ice  strength,  3) 


to  zero  ice  import  from  the  Arctic  Basin,  4)  to  zero 
currents,  S)  to  a  modified  current  field,  6)  to  zero 
winds  and,  finally,  7)  to  ice  dynamics  alone.  The 
results  of  these  simulations  are  discussed  in  the 
following  section. 


RESULTS  AND  DISCUSSION 

In  this  section  a  brief  discussion  of  the  current 
and  wind  fields  will  be  followed  by  the  various 
simulation  results.  Standard  simulation  results 
will  be  discussed  at  length,  followed  by  brief  dis¬ 
cussions  of  each  of  the  sensitivity  tests  which  com¬ 
pare  results  to  the  standard  run  and  to  observa¬ 
tions  where  possible,  in  the  Summary  and  Con¬ 
cluding  Remarks  section,  some  of  the  vital  results 
of  each  simulation  are  presented  in  tabular  form 
and  the  results  of  all  tests  are  summarized. 
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Figure  4.  Late  summer  currents,  average  ice  margins, 
and  manned  ice  station  tracks  (from  Einarsson  1972). 


Wind  and  current  fields 

The  60-day  averaged  wind  field  and  the  geostro- 
phic  current  field  for  the  October-November  1979 
period  of  study  are  shown  in  Figure  3.  The  most 
significant  feature  of  the  wind  field  is  the  narrow 
band  of  generally  northerly  winds  that  follows  the 
Greenland  coast.  The  surprising  elements  are  the 
large  topographical  influence  that  Greenland  ap¬ 
parently  has  on  the  surface  pressure  field,  and  the 
fact  that  this  feature  is  clearly  resolved  by  the 
NCC  analyzed  data  which  were  interpolated  from 
a  2 Vi -degree  latitude  and  longitude  grid.  As  a 
means  of  crude  verification,  the  pressure  fields 
were  manually  compared  to  those  produced  inde¬ 
pendently  by  Thorndike  and  Colony  (1980)  for  the 
same  period.  In  the  latter  analysis,  gridded  sea  le¬ 
vel  pressure  fields  were  constructed  using  an  opti¬ 
mal  interpolation  technique  applied  to  data  from 
approximately  13  drifting  buoys  in  the  Arctic  Ba¬ 
sin  and  70  high  latitude  land  stations.  The  manual 
comparison  showed  some  differences  between  the 
analyses,  namely  that  the  large  pressure  gradient 
along  the  Greenland  coast  was  not  as  apparent  in 
the  Thorndike  and  Colony  (1980)  analysis. 


Because  that  analysis  did  not  proceed  farther 
south  than  70  °N,  it  was  decided  to  use  the  NCC 
data. 

The  geostrophic  ocean  current  field  appears 
quite  smooth,  as  would  be  expected  from  the  fact 
that  a  temporally  constant  dynamic  height  field 
was  used  to  calculate  these  currents.  For  a  crude 
comparison,  late  summer  currents  as  compiled  by 
Einarsson  (1972)  are  shown  in  Figure  4.  The  geo¬ 
strophic  currents  have  the  same  general  direction 
as  those  of  Einarsson,  but  the  magnitudes  and  spe¬ 
cific  features  differ  considerably.  For  instance  the 
narrow  jet  of  high  velocity  currents  between  Spits¬ 
bergen  and  Greenland  is  shown  to  have  much 
higher  velocities  in  the  Einarsson  compilation. 
This  may  be  due  to  the  fact  that  some  of  Einars- 
son’s  data  represent  instantaneous  current  mea¬ 
surements  whereas  the  dynamic  heights  presum¬ 
ably  are  derived  from  a  long-term  data  base.  The 
question  remains,  however,  as  to  whether  geostro¬ 
phic  currents  are  representative  of  actual  currents, 
particularly  in  the  shallower  waters  adjacent  to  the 
Greenland  coast.  That  the  flow  here  may  be  par¬ 
tially  barotropic  rather  than  baroclinic  is  not  out 
of  the  question  (R.  Paquette,  pers.  comm.).  In  ad¬ 
dition,  it  is  well  known  that  the  motion  of  the  ice 
itself  transmits  stress  intw  the  ocean,  modifying 
the  currents  over  a  long  time  period.  These  prob¬ 
lems  concerning  actual  currents  can  only  be  re¬ 
solved  by  an  extensive  observation  network  or  a 
coupled  ice-ocean  model.  For  these  reasons,  the 
geostrophic  currents  which  seemingly  are  a  rea¬ 
sonable  first  order  approximation  of  the  currents 
in  this  area  are  used  in  this  study.  Simulations 
described  later  assess  the  response  of  the  model  to 
zero  currents  and  to  a  current  field  derived  from  a 
60-day  average  of  the  ice  velocity  field. 

Standard  simulation 

The  standard  simulation  represents  the  applica¬ 
tion  of  the  full  model  over  the  60-day  time  period, 
with  all  input  parameters  as  previously  described. 
The  idea  here  was  to  compare  the  model  results  to 
observations,  where  possible,  to  assess  the  overall 
validity  of  the  results.  In  addition,  these  “bench¬ 
mark”  results  are  used  for  comparison  to  other 
simulations  in  which  the  forcing  fields  are  varied . 

Initial  and  average  simulated  thickness  and 
compactness  fields  at  10-day  intervals  are  shown 
in  Figures  3  and  6.  The  ice  edge  positions  as  ob¬ 
tained  from  the  NPOC  ice  charts  for  times  closely 
corresponding  to  the  prediction  intervals  are  in¬ 
cluded  in  the  figures.  The  0.2  compactness  con¬ 
tour  (207*  concentration)  was  chosen  to  represent 


the  ice  edge  in  the  simulated  results.  Any  lower 
value  was  found  to  have  a  high  day-to-day  posi¬ 
tion  variability,  presumably  due  to  the  large  varia¬ 
tion  of  the  ice  growth  rates.  In  addition,  this  value 
appeared  to  correspond  well  with  the  0.1 -m  aver¬ 
age  thickness  contour,  and  both  seemed  to  be  rela¬ 
tively  stable  on  a  day-to-day  basis. 

These  figures  clearly  show  that  the  predicted  ice 
extent  is  excessive,  particularly  after  day  10.  Al¬ 
though  the  edge  as  indicated  on  the  NPOC  charts 
usually  enclosed  6-8  oktas  (concentration  in 
eighths),  the  predictions  are  still  excessive,  even  if 
a  higher  concentration  is  considered  as  the  pre¬ 
dicted  ice  edge  to  allow  for  possible  resolution  er¬ 
rors  when  the  NPOC  charts  were  compiled.  Some 
improvement  is  noted  on  day  60  if  the  predicted 
edge  is  taken  to  be  the  0.8  or  1 .0  compactness  con¬ 
tour,  however.  It  is  especially  evident  on  this  day 
because  the  predicted  compactness  is  much  more 
diffuse  than  on  previous  10-day  increments. 

A  more  quantitative  comparison  of  predicted 
and  actual  ice  extent  is  presented  in  Table  1.  In 
this  table  the  total  ice-covered  area  as  predicted  by 
the  model  is  compared  to  that  estimated  from  the 
ice  charts  for  the  10-day  intervals.  The  observed 
ice  coverage  was  determined  by  calculating  the 
product  of  the  area  covered  and  the  concentration 
specified  on  the  respective  ice  chart.  The  scale  of 
the  charts  plus  the  lack  of  detailed  compactnesses 
limits  the  accuracy  of  the  calculations;  however, 
the  overall  comparison  in  this  manner  is  felt  to  be 
meaningful.  Table  1  also  shows  the  percentage  dif¬ 
ference  in  predicted  versus  observed  coverage 
([predicted-observedj/observed)  and  the  percen¬ 
tage  change  for  both  predicted  and  observed  dur¬ 
ing  the  10-day  intervals. 

Table  1.  Amu  of  ice  cover  for  the  standard 

simulation  (10-day  intervals). 


Predicted 
(10"  m') 

Observed 
(10"  m‘) 

Difference 

(*) 

Initial  area 

1.80 

1.80 

Area  day  10 

1.88 

2.21 

-14.9 

Change  (4k) 

4.5 

22.7 

Area  day  20 

3.55 

2.24 

58.4 

Change  (4k) 

88.8 

1.4 

Area  day  30 

3.38 

2.06 

64.1 

Change  (4») 

-4.8 

-8.0 

Area  day  40 

3.79 

2.33 

62.7 

Change  (4k) 

12.1 

13.1 

Area  day  30 

4.49 

2.48 

81.0 

Change  (4k) 

18.5 

6.4 

Area  day  60 

5.84 

3.13 

80.8 

Change  (4k) 

30.1 

26.2 

Table  1  verifies  that  the  predicted  ice-covered 
areas  are  excessive  after  day  10,  and  the  simula¬ 
tion  ends  with  an  ice-covered  area  excess  of  81%. 
The  predicted  major  expansion  occurring  between 
days  10  and  20  is  primarily  due  to  large  growth 
rates  in  the  south  and  east,  and  this  will  be  exam¬ 
ined  further  in  subsequent  simulations.  It  is  inter¬ 
esting  to  note,  however,  that  after  this  period  the 
predicted  percentage  change  in  ice-covered  area 
tends  to  agree  with  that  of  the  observed.  Even  the 
decrease  in  ice  extent  between  days  20  and  30  is 
well  accounted  for  by  the  model. 

It  appears  that  the  high  growth  rates  are  primar¬ 
ily  responsible  for  the  large  ice  extent  that  is  pre¬ 
dicted.  This  reasoning  is  prompted  by  the  fact  that 
the  predicted  ice  edge  after  day  10  (Fig.  3  and  6)  is 
approximately  in  the  same  location  as  the  boun¬ 
dary  of  the  melt  rate  parameterization  discussed 
previously.  Once  the  ice  expands  to  this  limit,  then 
further  changes  in  extent  appear  to  be  due  to  a 
combination  of  dynamics  and  thermodynamics, 
but  the  magnitude  of  the  changes  is  limited  by  the 
ocean  heat  flux  parameterization,  at  least  until 
near  the  end  of  the  simulation  when  the  growth 
rates  are  high  enough  to  overcome  the  melt  rate 
specification.  The  upshot  is  that  the  growth  rates 
undoubtedly  need  to  include  a  better  parameter¬ 
ization  of  the  oceanic  heat  flux. 

Average  ice  velocities  are  also  useful  in  account¬ 
ing  for  thickness  and  compactness  variations. 
Ten-day  averaged  velocities  corresponding  to  the 
10-day  intervals  of  thickness  and  compactness  are 
shown  in  Figure  7.  Also,  the  60-day  averaged  vel¬ 
ocity  field  is  included  in  this  figure  (Fig.  7g). 
While  the  average  velocities  for  the  first  10  days 
are  nearly  negligible,  those  for  the  remainder  of 
the  10-day  intervals  closely  resemble  the  average 
wind  field  in  direction.  It  is  clear  that  ice  dynamics 
plays  some  role  in  the  large  ice  expansion  between 
days  1 1  and  20,  with  ice  being  advected  southward 
by  the  high  velocity  stream  (near  0.5  m  s  ')  shown 
in  Figure  7b.  As  the  prescribed  geostrophic  cur¬ 
rents  are  temporally  constant  and  no  larger  than 
0.05  m  s'1  (Fig.  3b),  this  large  velocity  increase 
over  the  first  10  days  can  be  attributed  to  winds. 
Although  temperature  fields  have  not  been  exam¬ 
ined  in  detail,  it  is  likely  that  these  winds  also  ad¬ 
vected  lower  air  temperatures  into  the  southern  re¬ 
gion  of  the  grid,  stimulating  a  rapid  ice  expansion. 
Ice  dynamics  may  also  have  been  partially  re¬ 
sponsible  for  the  decrease  in  ice  extent  between 
days  21  and  30.  The  average  velocities  for  this  pe¬ 
riod  (Fig.  7c)  show  a  marked  onshore  component 
in  the  vicinity  of  the  ice  edge.  This  velocity  config¬ 
uration  would  be  expected  to  confine  the  areal 


Figure  5  (corn’d).  Simulated  ice  thickness  fields. 
Dashed  line  is  observed  ice  edge  position  from 
NPOC  (1979). 


Figure  6  ( cont'd ).  Simulated  ice  compactness 
fields.  Dashed  line  is  observed  ice  edge  position 
from  NPOC  (1979). 


Figure  7.  A  veraged  simulated  ice  velocities. 


Figure  7  (cant ’d). 


Averaging  Interval  (days) 


56-60- - - r - : r- 

*  A  \  V 


t  t  ft  f- _ 


51-55 


V 


66-50,-- 

""STVT  v 


\  ^ 

'  V  ■  \  \  \  *  *  ^  *-»— » 

V 

*  ±  »  *  T 

v  ‘  « 


41-45*: 


♦  t  t 


36-40 


V 


■  0-  \  \  \  W\iM^ 


51-35* 


V  V  \  ^  ^  V 


26-30 


21-25 


V  V  i  v  i<  *  '  <  * 

V^TT'  ^  ^  '  * 

v-  ^V  ■  «■  -*  *  ■ 


\v^ 


1 

s' 

/ 

V  ,  ,  .  t  . 

- 

- 4 

*  V\ 

V 

IT 

v  * 

4  4 

4  4 

1-54  ♦  4- 

Coott 

-4  4  -  4  4  .4  1*  * 

40km 

t—  — 4 

as**-' 

X  x_ 

S#o  word 

Figure  8.  Five-day  averaged  velocities  for  the  row  of 
grid  points  immediately  above  the  11.0-m  coastal  ice 
buildup  shown  in  Figure  5g. 


coverage  by  advecting  ice  toward  the  coast.  In  ad¬ 
dition,  higher  air  temperatures  from  the  southeast 
could  be  expected  to  accompany  the  driving  wind 
field.  In  contrast,  the  average  velocities  for  days 
41  to  50  and  51  to  60  (Fig.  7e,  7f)  show  offshore 
velocity  components  in  this  region,  which  partially 
accounts  for  the  relatively  high  ice  expansion 
predicted  by  the  model  during  these  periods. 
Thus,  time  variations  in  the  wind  forcing  appear 
to  produce  time  variations  in  ice  extent,  both 
through  direct  forcing  by  ice  advection  and  more 
indirectly  through  variable  advections  of  air 
temperatures. 

The  ice  velocities  are  also  responsible  for  the 
predicted  thickness  and  compactness  variations 
occurring  near  shore  that  begin  to  become  obvious 
on  day  20  (Fig.  5c,  6c).  The  coastal  thickness 
buildups  and  areas  of  lesser  concentration  are  ob¬ 
viously  the  result  of  ice  dynamics  because  growth 
rates  are  very  small  for  ice  thicker  than  1 .0  m  and 
would  be  expected  to  have  little  effect  in  a  10-day 


period.  The  7.0-m  buildup  in  the  northern  section 
of  the  grid  on  day  20  (Fig.  5c)  continues  to  in¬ 
crease  to  1 1 .0  m  by  day  40  (Fig.  5e).  This  appears 
to  be  a  result  of  ice  impinging  upon  the  promon¬ 
tory  here  created  by  the  solid  boundary.  A  more 
detailed  look  at  the  ice  velocities  in  this  region  is 
presented  by  Figure  8,  which  shows  a  cross  section 
of  5-day  averaged  ice  velocities  for  the  row  of  grid 
points  immediately  above  this  large  buildup  and 
passing  through  the  entire  final  width  of  the  ice 
stream  in  this  region.  The  thickness  buildup,  re¬ 
sulting  from  ice  ridging,  occurs  just  beneath  the 
first  and  second  grid  points  from  the  left  in  this 
figure.  The  southeastward  drift  causes  the  ice  to 
accumulate  on  the  boundary  promontory  in  this 
region.  A  similar  figure  helps  explain  the  lesser 
concentration  area  which  occurs  immediately  to 
the  south  of  the  buildup  (Fig.  6c-g).  The  5-day 
velocity  averages  for  a  row  of  grid  points  beneath 
the  promontory  are  shown  in  Figure  9.  Once  again 
the  velocities  are  directed  southeastward.  In  this 
case,  however,  because  there  is  only  a  vertical 
boundary  at  the  coastline,  ice  is  advected  away 
from  the  coast,  eventually  resulting  in  lower  con¬ 
centrations  of  thinner  ice  (3.0  m). 

Figures  8  and  9  also  demonstrate  the  effects  of 
the  ice  strength  on  the  velocity  field.  In  Figure  8, 
the  coastal  buildups  have  increased  the  ice 
strength.  This  prevents  high  velocities  near  shore 
and  also  causes  the  ice  to  move  with  more  of  an 
offshore  component  into  areas  of  less  strength 
(smaller  thicknesses)  as  time  goes  on.  At  the  third 
grid  point,  strengths  are  much  less  due  to  lower 
thicknesses,  and  the  ice  moves  at  much  higher  vel¬ 
ocities  in  a  more  southerly  direction.  This  creates 
an  effective  velocity  shear  in  the  vicinity  of  the  ice 
buildup.  Further  offshore  the  velocity  shear  is 
more  likely  due  to  the  decreasing  winds  as  distance 
from  the  coast  increases.  In  contrast,  the  coastal 
velocity  shear  and  offshore  turning  are  barely 
discernible  in  Figure  9.  These  are  attributed  to  the 
fact  that  lower  strengths  are  maintained  adjacent 
to  the  coast  because  of  relatively  low  thicknesses 
and  compactnesses.  As  a  result,  velocities  are  not 
severely  affected. 

With  the  exception  of  the  excessive  ice  extent, 
the  predicted  thickness  and  compactness  fields  ap¬ 
pear  to  be  reasonable  and  consistent  with  the  spe¬ 
cified  input  fields  and  with  the  boundary  config¬ 
uration.  Whether  major  ridging  events  actually 
occur  in  the  predicted  thickness  buildup  locations 
is  unknown  because  there  are  no  data  available  for 
these  regions.  In  addition,  there  is  little  thickness 
information  available  for  this  entire  area.  Thick¬ 
nesses  in  the  northern  section  of  the  grid  generally 
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those  reported  by  Wadhams  (1980b), 
be  remembered  that  these  thicknesses 
lit  of  the  specified  initial  and  boundary 
Thicknesses  and  compactnesses  in- 
le  coast  is  approached  in  the  Denmark 
n  at  the  end  of  November  in  the  simu- 
Its  (Fig.  Sg  and  6g).  This  is  consistent 
s  obtained  from  an  analysis  of  sub- 
ar  data  collected  in  March  1971  in  this 
n  (Kozo  and  Diachok  1973,  Kozo  and 
4)  although  absolute  thicknesses  were 
d. 

uctive  to  examine  the  individual  roles 
port,  growth  and  export  have  on  the 
lume  in  this  standard  simulation.  Fig- 
ws  the  day-to-day  change  in  total  ice 
sng  with  the  volume  of  ice  produced 
r  growth  and  northern  inflow,  and  that 
it  of  the  southern  free  boundary.  Out- 
th  the  eastern  free  boundary  was  negli- 
figure  clearly  shows  that  thermody- 


Figure  10.  Simulated  daily  changes  in  total 
ice  volume  and  daily  volumes  of  growth, 
northern  inflow  and  southern  outflow. 


namics  (growth)  and  ice  dynamics  (northern  in¬ 
flow)  are  both  major  contributors  to  the  ice  mass 
balance  as  simulated  by  the  model.  Southern  out¬ 
flow,  as  expected,  only  contributes  to  the  balance 
during  the  latter  part  of  the  simulation  period.  Ice 
growth,  northern  inflow  and  southern  outflow 
correlate  with  the  daily  volume  change  with 
respective  coefficients  of  0.83,  0.85  and  0.30. 
These  coefficients  and  Figure  10  imply  that 
growth  and  inflow  had  nearly  equal  roles  in  the  ice 
mass  budget  according  to  the  simulation.  It  is  also 
interesting  to  note  that  a  0.48  correlation  coeffi¬ 
cient  exists  between  the  simulated  values  of  daily 
growth  and  northern  inflow.  Two  possible  explan¬ 
ations  may  account  for  this  correlation.  First,  dur¬ 
ing  periods  of  high  northern  inflow,  winds  would 
be  northerly,  advecting  colder  air  and  stimulating 
ice  growth  in  lower  concentration  areas.  Second, 
high  northern  inflow  is  likely  to  be  associated  with 
high  velocities  over  the  entire  grid.  This  would 
likely  create  areas  of  lesser  concentration  in  which 
new  ice  would  be  rapidly  produced.  The  large  vari¬ 
ability  in  northern  inflow,  including  the  reversals 
in  flow  direction,  indicate  that  the  simulated  ice 
transport  is  primarily  wind-induced. 
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Figure  ll.  Trajectories  of  ICEX  buoys  1564  and 
1568. 


A  note  of  caution  concerning  inflow  is  worth¬ 
while  at  this  point.  Inflow  is  partially  specified  as 
a  boundary  condition  by  the  fact  that  constant 
thicknesses  have  been  set  for  the  free  boundary 
cells  in  the  Fram  Strait.  By  way  of  comparison, 
Aagaard  and  Greisman  (1975)  estimated  that  the 
ice  outflow  rate  of  the  Arctic  Basin  was  approxi¬ 
mately  0.1  Sv  (3.154x10"  m’  yr'1).  The  total  in¬ 
flow  predicted  by  the  model  is  4.73x10"  m'  for 
the  two-month  simulation  period.  If  this  value 
were  to  remain  constant  for  the  entire  year,  2.86  x 
10"  m'  of  ice  would  be  imported  through  the 
Fram  Strait,  a  value  that  is  9.1%  less  than 
Aagaard’s  estimate.  The  simulated  outflow  may 
be  too  high,  however.  Einarsson  (1972),  in  an  in¬ 
vestigation  that  considered  only  area  of  ice  inflow 
(rather  than  volume),  calculated  that  73%  of  the 
total  yearly  ice  import  occurs  during  the  months 
of  December  through  May.  This  large  rate  was 
due  to  the  increased  northerly  winds  tl  at  typically 
occur  during  winter.  His  calculations  further 
showed  that  only  11%  of  the  total  inflow  took 
place  during  October  and  November.  If  this  is  in¬ 
deed  the  normal  case,  then  the  yearly  inflow  that 


the  model  would  predict  would  be  4.3  x  10"  m’. 
This  value  is  36%  higher  than  Aagaard’s  estimate. 
Because  these  *re  only  estimates  of  inflow,  how¬ 
ever,  it  is  difficult  to  assess  the  validity  of  the  sim¬ 
ulated  inflow.  The  simulated  value  at  least  seems 
to  be  of  reasonable  magnitude  for  the  two-month 
period. 

The  simulation  predicts  the  total  ice  volume  in¬ 
crease  to  be  7.85  x  10"  m’  for  the  two-month  peri¬ 
od.  This  constitutes  a  99%  increase  over  the  initial 
amount  of  ice.  Of  this  total  volume,  49.6%  was 
contributed  by  inflow  (with  southern  outflow  re¬ 
moved)  and  50.4%  was  added  by  growth.  That  in¬ 
flow  and  growth  produced  nearly  equal  volumes 
of  ice  is  somewhat  surprising.  The  suspicion  that 
too  much  ice  growth  is  taking  place  is  somewhat 
confirmed  by  once  again  comparing  relative  re¬ 
sults  to  those  of  Einarsson  (1972).  In  attempting 
to  establish  an  annual  budget  for  the  region  be¬ 
tween  76 °N  and  the  Denmark  Strait,  the  region 
where  most  of  the  predicted  ice  growth  takes 
place,  he  estimates  that  ice  growth  is  approximate¬ 
ly  'A  of  the  inflow  at  76  °N  for  the  months  of  Oc¬ 
tober  and  November.  In  this  region,  the  growth- 
to-inflow  ratio  predicted  by  the  model  is  at  least  as 
high  as  that  for  the  entire  grid  area.  The  predicted 
ratio,  which  is  1:1,  combined  with  an  ice  inflow 
rate  that  appears  reasonable,  indicates  that  far 
more  growth  is  taking  place  than  Einarsson  esti¬ 
mated.  Once  again,  it  appears  that  the  lack  of  a 
proper  parameterization  of  oceanic  heat  flux  is  al¬ 
lowing  excessive  ice  growth. 

The  accuracy  of  the  predicted  ice  velocities  at 
particular  locations  can  be  assessed  by  comparing 
them  to  the  velocities  of  buoys  that  were  drifting 
on  ice  floes  during  this  time  period.  The  trajecto¬ 
ries  of  ICEX  buoys  1564  and  1568  (Kloster  and 
Rafto  1980)  are  shown  in  Figure  1 1 .  An  interesting 
feature  of  the  trajectories  is  that  buoy  1568,  which 
is  located  some  50  to  140  km  closer  to  the  ice  edge, 
showed  a  much  larger  displacement  in  an  overall 
shorter  time  period  (8  days  less)  than  buoy  1564.  It 
is  also  apparent  that  the  speed  of  both  buoys  in¬ 
creases  southward.  Vinje  (1972,  1973,  1981)  and 
Wad  hams  (1981)  have  previously  reported  on  an 
acceleration  of  ice  as  distance  southward  and  dis¬ 
tance  from  the  coast  increase. 

To  compare  predicted  velocities  with  those  of 
the  buoys,  daily  velocities  were  interpolated  from 
the  grid  for  the  appropriate  buoy  location.  The 
predicted  and  observed  u  and  v  components  of 
velocity  are  shown  in  Figures  12  and  13.  Both  pre¬ 
dicted  and  observed  velocities  show  high  fre¬ 
quency  components  which  can  be  attributed  to  the 
fluctuating  winds.  It  appears  that  the  predictions 
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Figure  13.  Observed  and  predicted  components  of  velocity  for  K't'X  buoy  1568. 


Figure  14.  Actual  anti  predicted  trajectories  for  buoys  1564  and  1568.  Crosses  indicate  10-Julian-day  intervals 
(i.e.  300.  310.  ...). 


for  buoy  1564  arc  superior  to  those  for  buoy  1568. 
Additionally,  the  v  components  of  both  buoys 
seem  better  predicted  than  the  u  components.  Cor¬ 
relation  coefficients  between  predicted  and  ob¬ 
served  u  and  v  components  are  0.48  and  0.57,  re¬ 
spectively,  for  buoy  1564,  and  0.36  and  0.56  for 
buoy  1568.  Because  the  means  are  removed  when 
calculating  correlation  coefficients,  this  is  an  ef¬ 
fective  test  of  the  ability  of  the  model  to  predict 
only  the  high  frequency  components  of  the  buoy 
velocities.  This  is  demonstrated  by  the  high  corre¬ 
lation  coefficient  obtained  for  the  v  component  of 
buoy  1568,  where,  in  fact,  large  differences  be¬ 
tween  the  velocity  magnitudes  are  apparent.  These 
differences  become  clear  when  comparing  predict¬ 
ed  and  observed  means  for  buoy  1568  velocities. 
The  predicted  u  and  v  velocity  means  are  0.01  and 
-0.04  m  s'1  respectively,  while  those  for  the  ac¬ 
tual  velocities  are  0.15  and  -0.30  m  s'1.  In  con¬ 
trast,  the  u  and  v  predicted  velocity  means  for 
buoy  1564  are  0.08  and  -0.17  ms'1  and  those  of 
the  observed  velocities  are  0.08  and  -0.19  ms  '. 
Another  useful  statistic  for  assessing  the  predicta¬ 
bility  of  the  velocities  is  the  Root  Mean  Square  er¬ 


ror  (RMS  error  =  ['/n  I(Predicted-Observed)2] ,/l) 
wh'ch  gives  some  feeling  for  the  error  of 
amplitude  for  an  individual  velocity.  The  RMS  er¬ 
rors  for  the  u  and  v  velocities  were  0.12  and  0. 14  m 
s"1  for  buoy  1564  and  0.21  and  0.33  m  s"1  for  buoy 
1568.  Comparison  of  the  correlation  coefficients, 
velocity  means  and  RMS  errors  rapidly  verifies,  as 
do  Figures  12  and  13,  that  the  daily  velocities  of 
buoy  1564  are  more  accurately  simulated  than 
those  of  buoy  1568. 

The  accuracy  of  the  simulated  velocity  field  can 
also  be  assessed  by  calculating  the  simulated  tra¬ 
jectories  of  the  two  buoys.  The  trajectory  is  com¬ 
puted  by  interpolating  a  predicted  velocity  for  the 
predicted  buoy  position  of  the  previous  time  step, 
using  the  initial  buoy  position  for  a  starting  point. 
The  trajectories  for  both  buoys  are  shown  in  Fig¬ 
ure  14.  The  simulated  trajectories  are  not  satisfy¬ 
ing,  particularly  for  buoy  1568.  In  this  case  the 
calculated  trajectory  placed  the  buoy  within  the 
boundary  region  where  velocities  are  zero.  The 
simulated  trajectory  for  buoy  1564  is  somewhat 
better  but  its  "miss  distance"  for  the  final  buoy 
position  (day  334)  is  quite  large. 
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b.  Simulating  the  starting  position  90  km  west 
a.  Beginning  on  day  297.  of  the  actual  starting  position. 

Figure  15.  Actual  and  predicted  trajectories  for  buoy  1568. 


Simulated  trajectories  are  quite  sensitive  to  the 
time  and  place  chosen  for  the  buoy’s  initial  point. 
To  demonstrate,  the  extremely  poor  case  of  buoy 
1568  can  be  somewhat  improved  by  starting  the 
buoy  at  a  later  time  and  by  choosing  a  different 
starting  location.  Figure  15a  shows  a  simulated 
trajectory  obtained  by  starting  the  buoy  4  days 
later,  while  in  Figure  15b  the  starting  day  is  the 
original  day  (293)  but  the  starting  location  has 
been  moved  approximately  90  km  west.  Although 
the  simulated  trajectories  are  still  not  acceptable, 
significant  improvement  has  taken  place  over  the 
original  trajectory  shown  in  Figure  14b.  The  point 
to  be  made  here  is  that  small  errors  in  the  pre¬ 
dicted  velocity  field  over  only  a  few  days,  or  un¬ 
certainty  in  the  actual  buoy  starting  location,  can 
result  in  a  totally  unrealistic  trajectory.  This  can 
easily  lead  to  the  belief  that  velocities  in  the  region 
of  the  buoy  are  unrealistic  for  the  total  trajectory 
period. 

While  the  above  method  calculated  an  "ideal” 
trajectory,  the  results  are  very  sensitive  to  small  er¬ 


rors  in  the  simulated  velocity  field  or  initial  buoy 
location  as  has  been  shown.  Another  method  of 
computing  a  trajectory,  which  is  not  quite  so  sen¬ 
sitive,  is  to  again  sum  the  predicted  daily  veloci¬ 
ties,  only  this  time  to  take  these  velocities  from  the 
actual  daily  position  of  the  buoy,  rather  than  from 
the  predicted  position.  In  essence,  this  method 
consists  of  summing  the  predicted  velocities 
shown  in  Figures  12  and  13.  This  technique  is 
quite  useful  for  examining  the  long-term  cumula¬ 
tive  effects  (and  errors)  of  the  predicted  velocities 
at  the  locations  of  the  buoys.  Figure  16  shows  the 
trajectories  for  both  buoys  calculated  in  this  man¬ 
ner.  Once  again  it  becomes  clear  that  the  velocities 
predicted  for  buoy  1564  were  far  superior  to  those 
of  buoy  1568.  What  is  made  particularly  obvious 
in  this  figure  is  that  the  v  velocity  components 
simulated  for  buoy  1568  are  significantly  in  error, 
as  can  also  be  clearly  seen  from  Figure  13. 

The  behavior  of  these  trajectories  can  be  ana¬ 
lyzed  in  detail  using  the  10-day  averaged  velocity 
fields  (Fig.  7a-f);  however,  the  60-day  averaged 
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figure  16.  Actual  and  cumulative  daily  predicted  trajectories  for  buoys  1564  and  1568. 


velocity  field  can  be  used  for  a  brief  summary.  For 
descriptive  purposes,  the  actual  buoy  trajectories 
are  shown  superimposed  on  the  simulated  60-day 
averaged  ice  velocity  field  in  Figure  17.  The  trajec¬ 
tory  of  buoy  1564  placed  it  nearer  the  simulated 
high  ice  velocity  stream  adjacent  to  the  Greenland 
coast,  which  allowed  its  predicted  trajectories  (by 
either  method)  to  advect  it  southward.  The  predic¬ 
tions  of  buoy  1568  were  less  fortunate.  Its  starting 
position  (and  subsequent  positions)  placed  it  to 
the  east  of  this  predicted  high  velocity  region.  It  is 
for  this  reason  that  the  exercise  which  moved  its 
starting  location  90  km  west  significantly  im¬ 
proved  the  trajectory. 

The  surprising  feature  of  the  predicted  trajec¬ 
tories  is  that  the  poorest  predictions  were  for  buoy 
1568,  which  in  actuality  showed  the  largest  dis¬ 
placement  and  far  higher  velocities  than  buoy 
1 564.  This  results  from  the  fact  that  the  predicted 
velocities  closer  to  the  ice  edge  were  poor.  One 
major  problem  here  appears  to  be  related  to  the 
wind  fields.  A  comparison  of  the  60-day  averaged 
winds  (Fig.  3a)  and  the  60-day  averaged  ice  veloci¬ 
ties  (Fig.  7g)  leads  to  the  belief  that  the  high  veloc¬ 
ity  ice  stream  is  highly  dominated  by  the  winds, 
whose  magnitudes  may  be  too  large  near  the  coast 
and  too  small  at  the  ice  edge.  As  a  result,  the  ice 
velocities  fall  off  too  rapidly  to  the  east,  and  the 


on  the  60-day  averaged  ice  velocity  field. 
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a.  Thickness  field. 


b.  Compactness  field. 


Figure  18.  Sixty-day  thickness  and  compactness  fields  for  the  thermodynamic  simulation.  Dashed  line  is 
observed  ice  edge  position  for  2  December  1979  (NPOC  1979). 


high  velocities  that  apparently  existed  near  the  ice 
margin  during  this  period  are  not  properly  ac¬ 
counted  for. 

Other  processes  may  also  have  been  responsible 
for  the  observed  acceleration  of  the  ice  near  the 
edge.  Mesoscale  oceanographic  or  meteorological 
phenomena  associated  with  the  ice  edge  could  con¬ 
ceivably  cause  density  currents  or  winds  that 
would  not  be  resolvable  in  synoptic  scale  data.  In 
particular,  mesoscale  winds  resulting  from  local 
baroclinicity  or  katabatic  effects  may  be  occurring 
here.  Also,  temporally  varying  oceanic  processes 
would  not  be  accounted  for  by  the  climatological 
dynamic  height  Held  used  to  calculate  the  geo- 
strophic  currents  that  were  used  in  this  study. 

Another  possible  scenario  which  may  cause  an 
acceleration  of  ice  near  the  edge  has  recently  been 
suggested  by  Roed  and  O’Brien  (1981).  In  this 
work,  an  analysis  of  a  dispersive  medium  is  car¬ 
ried  out  in  which  a  momentum  equation  included 
a  pressure  term.  This  pressure  term  is  directed  nor¬ 
mal  to  the  ice  edge,  and  presumably  could  result 
from  the  random  bumping  of  ice  floes.  After  geo- 
strophic  adjustment,  the  velocity  field  exhibits  a 
jet-like  structure  near  the  ice  edge.  A  phenomenon 


of  this  nature  may  be  the  solution  to  the  problem 
presented  here  and,  if  so,  could  be  accounted  for 
by  the  model  used  in  this  study  by  a  modification 
of  the  constitutive  law  in  the  ice  margin  region. 
The  first  inclination,  however,  is  to  carry  out  a  de¬ 
tailed  examination  of  the  synoptic  wind  and  cur¬ 
rent  fields  used  in  this  study  to  see  if  they  ade¬ 
quately  represent  the  actual  winds  or  currents  in 
this  region. 

Because  all  other  simulations  provided  equally 
as  poor  or  worse  predictions  than  the  standard  run 
for  buoy  1568,  no  further  comparisons  to  this 
buoy  will  be  made.  In  addition,  when  trajectories 
are  calculated,  only  the  cumulative  daily  trajec¬ 
tory  method  will  be  used  because  it  points  out  all 
essential  features  and  because  the  “ideal”  trajec¬ 
tories  behaved  quite  poorly  on  the  remaining  sim¬ 
ulations,  being  so  sensitive  to  the  ice  velocities  at 
the  initialization  time  and  position. 

Thermodynamic  simulation 

For  this  simulation,  only  the  growth  rates  calcu¬ 
lated  for  the  region  were  of  concern.  This  sensitiv¬ 
ity  test  was  carried  out  simply  by  setting  all  ice 
velocities  to  zero  in  the  numerical  code.  In  this 
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manner,  thickness,  concentration  and  ice  extent 
can  increase  (or  decrease)  only  according  to  the 
growth  rates  derived  from  the  thermodynamic 
code. 

Figure  18  presents  the  thickness  and  compact¬ 
ness  fields  as  predicted  by  the  model  at  the  end  of 
the  60-day  simulation  period.  The  ice  extent  has 
increased  significantly  over  the  initial  field  (Fig.  Sa 
and  6a).  The  expansion  has  taken  place  almost  en¬ 
tirely  due  to  growth  of  thin  ice,  however.  The 
thicker  ice  (a  1.0  m)  has  expanded  very  little,  a 
good  indication  of  the  significant  difference  in  the 
growth  rates  between  thick  ice  and  thin  ice/open 
wsjser. 

Comparison  of  the  thickness  and  compactness 
fields  of  the  thermodynamic  simulation  (Fig.  18) 
to  those  of  the  60-day  thickness  and  compactness 
predicted  by  the  standard  run  (Fig.  5g  and  6g) 
gives  some  perception  of  the  effects  of  ice  dynam¬ 
ics.  The  most  obvious  difference  is  that  much 
more  thick  ice  occurs  near  the  coast  in  the  com¬ 
plete  simulation,  a  result  of  dynamically  induced 
advection  and  subsequent  ridging  of  ice  of  all 
thicknesses.  The  effects  of  divergence  which  cre¬ 
ated  the  lower  concentration  areas,  both  along  the 
coast  and  elsewhere,  are  also  not  apparent  in  the 
compactness  field  produced  by  the  thermody¬ 
namic  simulation. 

The  position  of  the  ice  edge  relative  to  the  re¬ 
ported  edge  position,  also  shown  in  Figure  18, 
shows  little  improvement  over  that  predicted  by 
the  standard  run.  This  tends  to  confirm  suspicions 
expressed  earlier  that  the  ice  growth  rates  were  re¬ 
sponsible  for  the  excessive  ice  extent.  This  is  par¬ 
ticularly  applicable  to  the  southern  section  of  the 
grid  where  the  largest  expansion  has  taken  place. 
To  examine  the  areal  extent  of  ice  in  more  detail,  it 
is  useful  to  prepare  a  table  similar  to  Table  1. 
Table  2  presents  the  area  of  ice  coverage  at  10-day 
intervals  along  with  the  percentage  increase  during 
the  intervals.  In  addition,  the  percent  difference  in 
coverage  between  the  thermodynamic  simulation 
and  both  the  reported  coverage  and  that  generated 
by  the  standard  run  are  included. 

The  table  clearly  shows  that  differences  of  areal 
ice  coverage  between  the  thermodynamic  and 
standard  simulations  are  small.  In  some  cases 
(days  30  and  40),  it  appears  that  the  ice  dynamics 
of  the  standard  run  were  acting  to  restrain  the  ice 
extent.  At  all  times,  however,  the  differences  are 
small  enough  that  the  excessive  ice  extent  can  be 
attributed  almost  entirely  to  thermodynamic 
growth.  Both  the  predicted  percentage  changes 
during  the  10-day  intervals  and  the  percent  dif¬ 
ference  from  the  observed  extent  are  similar  to 


Table  2.  Areas  of  ice  cover  for  the  thermodynamic 
simulation  (10-day  intervals). 


Predicted 
(10"  m') 

ft  m 

from  observed 

ft  diff  from 
sld  run 

Initial  area 

1.80 

Area  day  10 

1.82 

-17.6 

-3.2 

Change  (ft) 

l.l 

Area  day  20 

3.51 

56.7 

-1.1 

Change  (ft) 

92.8 

Area  day  30 

3.56 

72.8 

5.3 

Change  (ft) 

1.4 

Area  day  40 

3.85 

65.2 

1.5 

Change  (ft) 

8.1 

Area  day  50 

4.26 

71.8 

-5.4 

Change  (ft) 

10.6 

Area  day  60 

5.56 

77.6 

-4.8 

Change  (ft) 

30.5 

those  of  the  standard  run. 

That  the  thermodynamics  are  dominating  the 
ice  extent  and  edge  location  in  these  simulations 
should  not  be  construed  to  imply  that  ice  dynam¬ 
ics  are  not  relevant  to  this  process.  A  judgment  of 
this  nature  would  be  premature  at  this  point.  The 
problem  at  hand,  as  has  been  mentioned  several 
times  previously,  is  that  the  present  growth  rates 
seem  to  be  excessive.  This,  in  turn,  so  dominates 
the  simulated  ice  extent  that  the  effective  role  of 
dynamics  in  determining  ice  extent  appears  to  be 
small.  The  actual  role  of  the  ice  dynamics  in  this 
process  will  not  be  properly  resolved  until  more 
realistic  ice  growth  rates  are  utilized. 

The  total  volume  of  ice  produced  in  the  thermo¬ 
dynamic  insulation  is  much  less  than  that  pro¬ 
duced  in  the  standard  run.  The  volume  change 
here  was  2.57  x  10"  m’  compared  to  7.85  x  10"  m* 
for  the  previous  simulation.  In  this  simulation, 
growth  accounted  for  all  the  volume  change,  while 
in  the  standard  run  it  accounted  for  approximately 
half  of  the  change,  that  being  3.95  x  10"  m\  The 
standard  run,  then,  produced  53%  more  ice  by 
growth  than  the  thermodynamic  simulation.  This 
increase  in  growth  due  to  the  dynamics  is  likely 
due  to  new  ice  growth  in  areas  of  dynamically  in¬ 
duced  ice  divergence.  The  salient  point  here  is  that 
the  ice  dynamics  actually  increases  ice  production 
by  thermodynamics.  Similar  results  were  found  by 
Hibler  (1979)  in  a  modeling  study  of  the  Arctic 
Basin.  This  result  could  be  very  significant  to 
studies  dealing  with  air-sea  heat  exchange  in  this 
region.  The  implication  is  that  such  studies  must 
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Figure  19.  Sixty-day  thickness  and  compactness  fields  for  the  zero-strength  simulation.  Dashed  line  is  observed 
ice  edge  position  for  2  December  1979  (NPOC  1979). 


include  the  effects  of  ice  dynamics  to  properly 
treat  air-sea  energy  exchanges. 

Zero  ice  strength 

The  effect  of  the  internal  ice  stress  term  in  the 
momentum  balance  can  be  assessed  by  allowing 
the  ice  to  have  no  strength.  This  damps  out  ice  in¬ 
teraction  with  itself  and  effectively  creates  a  free 
drift  situation.  With  the  zero  strength  condition 
imposed,  the  ice  also  has  no  resistance  to  deforma¬ 
tion.  In  practice,  this  case  is  simulated  by  setting 
P*,  an  empirical  constant  in  eq  4,  to  zero.  This  ef¬ 
fectively  results  in  the  ice  strength  and  bulk  and 
shear  viscosity  terms  all  being  zero,  thus  eliminat¬ 
ing  the  internal  ice  stress  term  £  in  eq  1 . 

Thickness  and  compactness  fields  at  the  end  of 
the  60-day  simulation  period  are  shown  in  Figure 
19.  A  salient  characteristic  of  the  zero  strength 
condition  is  manifested  by  the  unreasonably  large 
thickness  buildups  that  have  occurred  adjacent  to 
the  coast.  The  necessity  of  allowing  the  ice  to  in¬ 
teract  with  itself  in  any  effort  to  model  this  region 
is  clearly  demonstrated  by  this  figure.  Farther 
east,  nearer  the  ice  edge,  thicknesses  appear  to  be 


more  reasonable;  thus  it  seems  that  a  free  drift 
model  may  perform  adequately  here. 

The  60-day  averaged  velocity  field  for  the  zero 
strength  simulation,  shown  in  Figure  20,  helps  ex¬ 
plain  some  of  the  features  apparent  in  the  thick¬ 
ness  and  compactness  fields.  A  definite  onshore 
velocity  component  is  obvious  in  several  locations 
along  the  coast.  This  onshore  component,  which  is 
a  result  of  strength  being  independent  of  thickness 
(actually  zero  in  this  case),  amplifies  the  effect  cre¬ 
ated  by  the  ice  having  no  resistance  to  deforma¬ 
tion  and  results  in  the  physically  unrealistic  thick¬ 
nesses.  Comparison  of  this  velocity  field  to  that  of 
the  standard  run  (Fig.  7g)  manifests  the  rectifying 
effect  that  ice  interaction  allows.  Higher  ice  stress¬ 
es  near  the  coast  in  the  standard  run  effectively 
curtail  further  motion  in  that  direction,  yet  allow 
motion  toward  areas  of  less  stress. 

The  compactness  “holes”  which  appear  short 
distances  away  from  the  coast  are  also  created  by 
the  ice  dynamics.  Here,  ice  is  simply  being  ad- 
vecied  away  from  a  particular  area  faster  than  it 
can  be  replaced  by  advert  ion  from  an  adjacent 
area  or  by  growth.  This  results  in  the  creation  of 
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Figure  20.  Sixty-day  averaged  velocity  field  for  the 
zero-strength  simulation. 


low  concentration  cells.  The  phenomenon  is  pre¬ 
vented  in  the  standard  simulation  because  velocity 
magnitudes  near  the  coast  are  decreased  by  higher 
strengths.  With  these  numerous  low  concentration 
areas,  it  is  not  surprising  that  the  area  of  ice 
coverage  at  the  end  of  the  simulation  is  4.47  x  10' 1 
m\  a  value  that  is  23%  less  than  the  standard  run. 

The  total  ice  growth  for  this  simulation  was 
1.33  x  10"  m\  a  value  that  is  two  orders  of  magni¬ 
tude  higher  than  that  of  either  the  standard  run  or 
the  thermodynamic  simulations.  This  unreason¬ 
ably  high  value  is  presumably  due  to  the  growth  of 
ice  in  the  areas  of  ice  divergence.  The  reasoning 
here  is  that  if  nearly  all  the  ice  is  advected  out  of 
these  areas  (the  low  compactness  cells)  at  every 
time  step,  the  high  growth  rates  of  thin  ice  and 
open  water  will  be  continually  sustained,  leading 
to  excessively  high  total  ice  growth.  In  contrast, 
the  total  inflow  for  the  period  appears  to  be  more 
reasonable,  being  6.18x10"  m'.  This  is  30% 
larger  than  the  inflow  predicted  by  the  standard 
run.  This  is  not  surprising,  however,  because  high¬ 
er  velocities  arc  to  be  expected  for  the  inflow  re- 


Figure  21.  Actual  and  cumulative  daily  predicted 
trajectories  for  buoy  1564  for  the  zero-strength  si¬ 
mulation. 


gion,  and  these  are  indeed  in  evidence  in  Figure 
20. 

The  correlation  coefficients  between  the  pre¬ 
dicted  and  observed  velocities  of  buoy  1 564  were 
0.45  for  both  u  and  v  components.  These  coeffi¬ 
cients  are  less  than  those  of  the  standard  run,  par¬ 
ticularly  with  the  v  component  (0.57  for  the  stan¬ 
dard  run).  The  predicted  velocity  means  are  -0.03 
and  -0.17  m  s"  for  the  u  and  v  components  (0.08 
and  -0.19  m  s"1  for  the  observed).  The  u  compo¬ 
nent  mean  shows  an  excessive  onshore  velocity 
trend.  RMS  errors  are  also  higher  than  those  of 
the  standard  run,  being  0.16  m  s'1  for  both  u  and  v 
components  (0.12  and  0.14  m  s"  for  the  standard 
run).  The  true  deficiencies  in  the  velocity  predic¬ 
tions,  however,  are  made  more  apparent  by  the 
cumulative  daily  predicted  trajectory  for  buoy 
1564,  which  is  shown  in  Figure  21.  The  excessive 
onshore  component  of  velocity  is  quite  evident  in 
this  figure,  eventually  placing  the  buoy  well  into 
the  boundary  area  representing  Greenland.  Once 
again,  the  lack  of  a  velocity  rectification  effect 
produced  by  ice  interaction  is  made  clear  by  the 
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Figure  22.  Sixty-day  thickness  and  compactness  fields  for  the  zero-inflow  simulation.  Dashed  line  is  observed 
ice  edge  position  for  2  December  1979  (NPOC  1979). 
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comparison  of  this  trajectory  with  that  of  the  stan- 
dard  run  (Fig.  16a). 

Zero  ice  import 

Previous  simulations  have  shown  that  the  ice  in¬ 
flow  through  the  Fram  Strait  constitutes  a  major 
part  of  the  mass  budget  in  this  region.  In  this  light, 
a  worthwhile  sensitivity  test  is  to  not  allow  inflow 
and  assess  the  impact  of  this  on  the  model  results. 
Zero  inflow  is  simulated  simply  by  specifying  zero 
thicknesses  for  the  northern  free  boundary  cells. 

Figure  22  shows  the  thickness  and  compactness 
fields  at  the  end  of  the  60-day  simulation.  These 
fields  are  similar  to  the  analogous  fields  of  the 
standard  run  (Fig.  5g  and  6g),  with  several  notable 
exceptions.  These  discrepancies  are  particularly 
apparent  in  the  northern  section  of  the  grid.  The 
most  apparent  difference  is  that  thicknesses  are 
generally  lower  in  this  section.  Coastal  buildups  in 
this  region  are  much  less  than  those  produced  in 
the  standard  simulation.  Similar  features  occur  in 
the  compactness  field  where  lower  concentrations 
have  developed  in  the  northern  section.  The  south¬ 
ern  third  of  the  grid  area  and  the  total  ice  extent 


are  virtually  unchanged  between  the  two  simula¬ 
tions,  however.  The  cause  for  these  differences  in 
the  northern  sector  is  obviously  that  ice  is  being 
transported  southward  or  ridged  near  the  coast, 
and  no  thicker  ice  is  being  advected  in  through  the 
free  boundary  to  replace  it.  Also,  new  ice  growth 
is  not  sufficient  to  sustain  the  100%  concentration 
level  in  this  high  velocity  region.  However,  the 
growth  of  thin  ice  is  occurring  at  a  relatively  high 
rate  due  to  the  advection  of  ice  out  of  the  northern 
region.  The  total  ice  growth  for  this  simulation 
was  5.62  x  10"  m',  a  40%  increase  over  the  stan¬ 
dard  run.  Because  velocities  are  less  near  the  ice 
edge,  its  position  is  predominantly  controlled  by 
thermodynamics,  as  was  the  case  for  the  standard 
simulation. 

That  the  southern  region  of  the  simulation  area 
is  free  from  the  effects  of  no  inflow  during  this 
simulation  is  verified  by  the  volume  of  ice  exiting 
the  southern  open  boundary.  For  both  the  stan¬ 
dard  run  and  the  zero  inflow  cases,  the  total 
southern  outflow  volume  for  the  60  days  is  0.83 
x  10"  m'.  The  implication  is  that  during  this  peri¬ 
od,  the  effects  of  ‘‘running  out"  of  thick  ice,  as 
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Figure  23.  Sixty-day  averaged  velocity  field  for  the 
zero-inflow  simulation. 


has  occurred  in  the  north,  had  not  yet  reached  the 
southern  outflow  region.  This  fact  is  obvious  from 
the  similarity  of  the  thickness  and  compactness 
fields  to  those  of  the  standard  run  in  the  southern 
region.  The  indication  is  that  the  southern  region 
would  eventually  be  affected,  but  an  appreciably 
longer  simulation  would  be  required  to  sort  out 
these  effects. 

The  60-day  averaged  velocity  field  is  presented 
in  Figure  23.  With  the  exception  of  the  northern 
region,  where  velocities  have  a  slightly  greater  on¬ 
shore  component,  the  velocities  are  nearly  identi¬ 
cal  with  those  of  the  standard  run  (Fig.  7g).  The 
larger  onshore  component  and  slightly  higher 
magnitudes  can  be  attributed  to  lower  ice 
strengths  in  the  region.  Here  lower  strengths  have 
again  been  induced  by  the  lower  ice  thicknesses 
and  concentrations. 

Correlation  coefficients  of  the  u  and  v  velocity 
components  with  those  of  buoy  1364  are  0.48  and 
0.52,  respectively.  The  v  coefficient  is  only  slightly 
less  than  that  of  the  standard  run  (0.57).  The  pre¬ 
dicted  velocity  means  are  0.06  and  -0.16  m  s''  for 
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Figure  24.  Actual  and  cumulative  daily  predicted 
trajectories  for  buoy  1564  for  the  zero-inflow  si¬ 
mulation. 


the  u  and  v  components— only  slightly  different 
from  those  of  the  standard  run  (0.08  and  -0. 19  m 
s"').  RMS  errors  for  the  predicted  components  are 
0.14  and  0.15  m  s'1,  once  again  close  to  those  of 
the  standard  run  (0.12  and  0.14  m  s  ").  The 
cumulative  daily  predicted  trajectory  for  this  buoy 
is  shown  in  Figure  24.  The  slight  onshore  velocity 
component  is  again  emphasized  in  this  figure,  with 
the  buoy  being  placed  slightly  closer  to  shore  than 
with  the  standard  simulation.  The  trajectory  in 
this  case  also  misses  the  final  buoy  position  by  a 
larger  distance  than  the  standard  run,  bearing  out 
the  lower  correlation  coefficient  and  less  negative 
mean  for  the  v  component.  Overall,  however,  the 
velocities  do  not  seem  greatly  affected  by  the  zero 
inflow  stipulation.  A  longer  simulation,  in  which 
thicknesses  decreased  over  the  entire  region,  might 
show  significant  effects. 

Zero  currents 

This  sensitivity  test  examines  the  effect  of  the 
geostrophic  currents  specified  for  the  other  simu¬ 
lations  simply  by  turning  those  currents  off.  With 
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Figure  25.  Sixty-day  averaged  ice  velocities  for  the 
zero-current  simulation. 


(his  specification,  (he  ice  can  be  (hougtu  of  as 
moving  across  a  stagnant  ocean.  Water  stress  is 
still  an  integral  part  of  the  momentum  equation, 
only  it  is  calculated  with  a  zero  current  velocity. 
The  force  due  to  the  tilt  of  the  sea  surface  is  also 
zero  in  this  case. 

At  the  end  of  the  60-day  simulation  period,  the 
thickness  and  compactness  fields  were  nearly  iden¬ 
tical  to  those  of  the  standard  run;  thus  they  are  not 
shown  here.  Even  regions  of  coastal  buildups  and 
low  concentrations  had  approximately  the  same 
thickness  and  compactness  values.  This  is  indica¬ 
tive  of  the  nearly  negligible  effect  exhibited  by  the 
geostrophic  currents  specified  for  the  previous 
simulations  on  the  ice  dynamics.  The  inference  is 
that  the  ice  dynamics  in  previous  simulations  has 
been  primarily  wind-driven. 

With  such  obvious  similarities  in  the  thickness 
and  compactness  fields,  it  is  not  surprising  that  the 
total  area  of  ice  coverage  at  the  end  of  60  days  is 
5.78  x  10"  mJ.  This  is  within  10%  of  the  value  pre¬ 
dicted  by  the  standard  run.  The  total  volume  of  ice 
exiting  the  southern  boundary  is  also  quite  similar 


Figure  26.  Actual  and  cumulative  daily  predicted 
trajectories  for  buoy  1564 for  the  zero-currents  si¬ 
mulation. 


to  that  of  the  standard  run,  being  0.81  x  10"  m' 
(0  83x10"  m'  for  the  standard  run).  Differences 
are  apparent  in  the  volume  of  northern  inflow, 
however.  In  the  zero  current  simulation  that 
volume  is  3.85x10"  m\  approximately  20%  less 
than  the  inflow  of  the  standard  simulation. 

The  differences  in  the  volume  of  northern  in¬ 
flow  ice  can  be  attributed  to  the  relatively  high 
current  velocities  in  th?  inflow  area,  as  shown  in 
Figure  3b.  This  region  would  then  be  expected  to 
be  more  severely  influenced  when  currents  are  set 
to  zero.  This  is  manifested  by  the  60-day  averaged 
ice  velocities,  shown  in  Figure  25.  Here,  velocities 
of  ice  entering  the  grid  are  shown  to  be  slightly  less 
than  those  of  the  analogous  figure  for  the  stan¬ 
dard  run  (Fig.  7g).  In  addition,  this  field,  which  is 
essentially  the  result  of  dynamics  driven  primarily 
by  wind,  shows  a  small  velocity  component  that 
would  drive  ice  out  of  this  boundary,  thus  further 
reducing  the  net  inflow  volume.  Further  south,  the 
velocity  vectors  appear  to  be  identical  to  those  of 
the  standard  simulation. 

In  further  assessment  of  predicted  ice  velocities. 
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the  correlation  coefficients  of  the  u  and  v  velocity 
components  with  those  of  buoy  1564  are  0.47  and 
0.54.  As  with  other  simulations,  these  are  quite 
similar  to  those  of  the  standard  run,  only  being 
slightly  smaller  in  the  v  component.  The  RMS  er¬ 
rors  are  exactly  the  same  as  those  of  the  standard 
run,  being  0.12  and  0. 14  ms1  for  (he  u  and  v  com¬ 
ponents.  The  mean  u  and  v  velocity  components, 
0.07  and  -0.15  m  s'1,  are  again  comparable  to  the 
observed  means  (0.08  and  -0.19  m  s'1),  and  quite 
close  to  those  of  the  standard  run  (0.08  and  -0.17 
m  s'1).  The  cumulative  daily  predicted  buoy  trajec¬ 
tory  is  shown  in  Figure  26.  The  predicted  trajec¬ 
tory  is  slightly  further  west  than  that  predicted  by 
the  standard  run,  probably  due  to  the  greater  on¬ 
shore  component  of  velocities  near  the  southern 
inflow  region. 

In  general,  the  geostrophic  currents  specified 
for  the  other  simulations  had  little  effect  on  the 
model  results  except  at  the  inflow  region.  The 
most  significant  effect  was  to  increase  the  volume 
of  ice  entering  the  region  by  contributing  to  larger 
southward  velocities.  Thicknesses  and  concentra¬ 
tions  over  the  remainder  of  the  grid  seem  unaffect¬ 
ed  by  turning  off  the  geostrophic  currents. 

Modified  currents 

Because  the  previous  simulation  showed  that 
the  geostrophic  currents  had  so  little  impact  on  the 
model  results,  it  was  decided  to  dramatically  alter 
the  current  field.  This  step  was  taken  partially  be¬ 
cause  previous  investigators  have  attributed  a  ma¬ 
jor  component  of  the  ice  transport  in  this  region  to 
currents  (Vowinckel  1964,  Einarsson  1972).  This 
study  is  somewhat  limited,  however,  because  the 
effects  of  current  transport  can  only  be  examined 
for  the  60-day  simulation  period.  In  another 
modeling  study  of  the  arctic  ice  cover  far  from 
shore  (Hibler  and  Tucker  1979),  currents  had  a 
negligible  short-term  effect  on  ice  drift  but  were 
found  to  be  important  in  the  long-term  drift.  As  a 
result,  the  feeling  was  that  only  a  larger  magnitude 
current  field  would  show  a  significant  effect  on 
the  ice  dynamics  for  this  relatively  short-term 
study. 

For  this  simulation,  the  60-day  averaged  ice  vel¬ 
ocity  field  generated  by  the  zero  current  simula¬ 
tion  (Fig.  25)  was  used  as  a  temporally  constant 
current  field.  The  basic  idea  behind  this  was  to 
simulate  a  barotropic  oceanic  flow  in  the  East 
Greenland  area,  where  currents  would  be  a  two- 
month  average  of  the  ice  velocities.  This  situation 
is  probably  not  realistic  because  there  is  only  a  rel¬ 
atively  narrow  region  of  shallow  water  adjacent  to 
the  coast,  and  other  topographic  features  in  the 


area  (subsea  ridges  and  sills)  would  not  be  amen¬ 
able  to  barotropic  flow.  In  addition,  it  is  not  clear 
that  ocean  currents  beneath  ice  covers  are  entirely 
driven  by  stress  transmitted  into  the  ocean  from 
the  moving  ice.  Another  reason  that  the  averaged 
ice  velocity  field  was  used  to  simulate  the  steady 
current  field  was  because  the  velocities  are  gener¬ 
ally  an  order  of  magnitude  higher  than  those  of 
the  geostrophic  currents,  and  the  sensitivity  of  the 
model  to  much  larger  currents  was  of  interest. 

The  60-day  averaged  ice  velocity  field  is  shown 
in  Figure  27.  The  effect  of  the  increased  current 
velocities  is  immediately  apparent.  Here,  ice  veloc¬ 
ities  are  much  higher  than  (hose  of  the  standard 
run.  In  some  cases,  particularly  in  the  area  adja¬ 
cent  to  the  coast,  velocities  arc  two  to  three  times 
larger  than  those  predicted  by  the  standard  run 
(Fig.  7g).  There  appears  to  be  no  major  direction 
change,  however,  except  in  eastern  portions  of  the 
northern  inflow  region  where  the  velocities  exhibit 
a  much  larger  onshore  component.  Because  this 
was  a  characteristic  of  the  current  field  used  here, 
this  is  not  surprising. 

The  60-day  thickness  and  compactness  contours 


Figure  27.  Sixtv-day  averaged  velocity  field  for  the 
modified-currents  simulation. 
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av-day  thickness  and  compactness  fields  for  the  ntodified-currents  simulation.  Dashed  line  is 
•dye  position  for  2  December  1979  (NPOC  1979). 


in  Figure  28.  Only  small  differences 
lese  fields  from  those  produced  by 
run  (Fig.  5g  and  6g).  The  coastal 
ch  occur  in  the  same  locations  as  the 
are  slightly  larger.  In  addition,  areas 
enlration  evident  in  the  standard  run 
er  compactnesses  in  this  simulation, 
h  be  attributed  to  (he  higher  ice  vel- 
ncreased  drift  can  be  expected  to 
igher  thicknesses  on  boundary  pro- 
ihe  drift  path  and  to  lower  ice  con- 
their  lee,  where  velocities  arc  mov- 
rom  (he  coast. 

Hive  daily  predicted  trajectory  lor 
hown  in  Figure  29.  The  effect  of  the 
H'ities  is,  again,  obvious  in  this  fig- 
predicted  final  position  of  the  buoy 
st  of  the  actual  position.  The  pre- 
v  components  of  velocity  correlate 
the  buoy  with  coefficients  of  0.42 
lectively.  These  are  somewhat  less 
the  standard  simulation  (0.48  and 
nedly  due  to  the  higher  velocities, 
leans  significantly  reflect  these  Itigh- 
ieing  0.13  and  -0.25  ms1  for  the  it 


and  v  components  (0.08  and  -0.17  m  s  '  for  the 
standard  run).  This  is  the  only  simulation  in  which 
the  magnitudes  of  the  velocity  means  are  larger 
than  the  observed  component  means  (0.08  and 
-0.19  m  s  ').  Likewise,  the  RMS  errors  between 
the  predicted  and  observed  velocity  components 
are  larger  than  in  any  previous  simulation  (0.20 
and  0.21  m  s  ’). 

The  total  area  of  ice  cover  at  the  end  of  the  sim¬ 
ulation  period,  5.88x10"  m2,  is  very  similar  to 
that  of  the  standard  simulation  (5.84x10"  m2). 
This  is  not  surprising  because,  as  previously  noted, 
thermodynamics  seems  to  dominate  the  areal  ice 
coverage  unless  significant  divergence  due  to  ice 
dynamics  is  taking  place.  The  inflow  and  outflow 
volumes  do  show  large  differences,  however.  The 
northern  inflow  of  7.41x10"  m'  is  50%  higher 
than  that  of  the  standard  run.  The  southern  out¬ 
flow,  2.94x10"  m\  is  more  than  three  times  that 
of  the  standard  run.  In  addition,  the  total  growth 
for  the  period  was  5.61  x  10"  m\  some  40%  higher 
than  the  growth  that  occurred  in  the  standard  sim¬ 
ulation.  This  increase  in  growth  may  also  be  at¬ 
tributed  to  the  increased  ice  velocities,  which,  in 
stimulating  the  overall  ice  dynamics,  created  more 
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Figure  29.  Actual  and  cumulative  daily  predicted 
trajectories  for  buoy  1564  for  the  modified-cur- 
rents  simulation. 


areas  of  divergence  in  which  new  growth  took 
place. 

These  results  show  that  the  model  is  indeed  very 
sensitive  to  the  current  field  in  this  region.  The 
predicted  increase  in  overall  ice  transport  is  rather 
expected  from  the  order  of  magnitude  increases  in 
current  velocities.  The  results  tend  to  indicate, 
however,  that  the  current  velocities  used  in  this 
simulation  are  probably  too  large.  This  is  implied 
primarily  by  the  predicted  buoy  drift  trajectory, 
which  places  the  buoy  too  far  to  the  southeast, 
and  by  the  excessive  inflow  and  outflow  volumes. 
Of  course,  this  line  of  thought  assumes  that  all 
other  model  parameters,  including  the  air  and 
water  stress  drag  laws,  are  reasonable.  An  appro¬ 
priate  current  field  for  this  region  will  presumably 
not  be  available  until  coupled  ice-ocean  model 
studies  are  undertaken .  What  would  appear  to  be 
a  more  appropriate  current  field  at  this  point  is 
something  between  the  geostrophic  field  used  in 
previous  simulations  and  the  averaged  ice  velocity 
field  used  here. 


Zero  winds 

The  idea  of  this  simulation  was  to  assess  the  rel¬ 
ative  importance  of  the  winds  as  a  driving  force  in 
this  short-term  study.  It  has  been  previously  con¬ 
cluded  that  the  geostrophic  currents  used  in  other 
simulations  have  only  a  small  impact  on  the  thick¬ 
ness  and  velocity  fields.  In  addition,  a  test  run  in 
which  winds  were  set  to  zero  and  the  geostrophic 
currents  alone  forced  the  ice  dynamics  showed  re¬ 
sults  quite  similar  to  the  thermodynamics  simula¬ 
tion.  Ice  velocities  were  very  small,  resulting  in  a 
northern  inflow  volume  that  was  only  10%  that  of 
the  standard  simulation.  Southern  outflow  was 
nonexistent.  Because  it  is  suspected  that  the 
geostrophic  currents  may  not  be  representative  of 
actual  currents,  the  modified  current  field  of  the 
previous  simulation  in  which  the  current  field  is 
the  60-day  averaged  ice  velocity  field  generated  by 
the  zero  current  simulation  has  been  used  here. 
The  current  velocities,  as  previously  noted,  are 
suspected  to  be  somewhat  excessive.  In  this  light, 
this  simulation  provides  a  sort  of  extreme  test  of 
the  influence  of  temporally  constant  currents  on 
the  model  results. 

Figure  30  presents  the  60-day  thickness  and 
compactness  fields  for  this  zero  wind,  modified 
currents  simulation.  Here,  minor  effects  of  the 
currents  on  ice  dynamics  are  apparent.  In  particu¬ 
lar,  the  3.0-m  contour  at  the  inflow  boundary  and 
the  slightly  lower  concentrations  along  the  coast 
are  current  effects.  Also,  the  ice  cover  here  is 
slightly  more  expansive  than  that  generated  by  the 
standard  or  thermodynamics  simulations.  The 
total  areal  coverage  is  6. 15  x  10' '  m\  a  5%  increase 
over  the  standard  run  and  10%  greater  than  the 
thermodynamics  simulation.  That  the  ice  veloci¬ 
ties  stimulate  this  slight  expansion  is  apparent 
from  the  60-day  averaged  velocity  field  shown  in 
Figure  31.  The  expansion  is  largely  in  the  northern 
sector  of  the  grid,  and  it  is  here  that  velocities  have 
more  of  an  easterly  component  than  previous  sim¬ 
ulations  have  shown. 

The  velocity  field,  in  particular  the  stream  adja¬ 
cent  to  the  coast,  shows  magnitudes  that  are  near¬ 
ly  one-half  those  of  the  standard  simulation.  This 
is  also  reflected  in  the  volumes  of  inflow  and  out¬ 
flow  ice,  wviose  values  are  2.02x10"  m'  and 
0.45x10"  m\  respectively.  These  are  approxi¬ 
mately  half  those  of  the  standard  run,  which  used 
the  lower  velocity  geostrophic  currents.  These  val¬ 
ues  are  also  about  50%  of  those  produced  by  the 
zero  current  simulation.  This  indicates  that  during 
this  period,  the  model  predicts  that  ice  transport 
by  currents  (which  are  fell  to  be  excessive)  is  on 
the  order  of  half  of  that  transported  by  wind. 
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Figure  30.  Sixty-day  thickness  and  compactness  fields  for  the  zero-winds  simulation.  Dashed  tine  is  observed 
ice  edge  position  for  2  December  1979  (NPOC  1979). 

Ji  l  l  Expectedly,  the  cumulative  daily  predicted  drift 

r-T*  *  *  trajectory,  shown  in  Figure  32,  leaves  buoy  1564 

— J"!  '  l  .  far  shy  of  its  final  observed  position.  This  is  veri- 

. .  fied  by  the  v  velocity  component  mean,  -0.07  m 

■» » •  •  s'1,  which  has  the  lowest  magnitude  of  any  simula- 

»*.  tion  and  indicates  far  too  little  southward  trans- 

I  *  *  port  when  compared  to  the  observed  v  component 

|  ■  ■  mean  (-0.19  m  s*’).  The  predicted  u  component 

•  mean  (0.06  m  s'1),  on  the  other  hand,  is  quite  com- 

,  parable  to  the  observed  mean  (0.08  m  s'1).  Surpris- 

vwv* >  ’  •  •  •  •  ingly,  the  u  component  RMS  error  is  0.09  ms'1, 

*  lower  than  any  previous  simulation,  while  the  v 

v  •  □  component  RMS  error,  0.17  m  s"',  is  somewhat 

large.  Comparison  of  the  predicted  and  observed 

{  ?  ^  i  *  * . ♦  u  component  velocities  (not  shown  here)  revealed 

i£\\'*’*'*****  that  Pre£*‘cteci  velocity  was  nearly  constant  due 

LAnn ••••***  to  the  lack  of  fluctuations  caused  by  the  winds. 

Yssw'  ******  Because  the  observed  magnitudes  of  the  buoy  vel- 

\  \  ocity  u  component  are  relatively  small  (Fig.  12), 

{*/*::::  the  cumulative  squared  differences  used  in  the  er- 

N  <  *  *  *  r  ror  calculation  were  less  than  for  other  simula- 

\jl***4*  tions,  resulting  in  a  smaller  RMS  error  value.  In 

J  contrast,  the  correlation  coefficients  between  the 

s  *  '*  i  predicted  u  and  v  velocity  components  and  those 

C  of  the  buoy  are  0.03  and  0.29.  These  coefficients, 

— I.  J  which  are  the  lowest  of  any  of  the  simulations 

which  included  ice  dynamics,  emphasize  the  day- 
Figure3l.  Sixty-day  averaged  velocity  field  for  the  to-day  variation  in  buoy  velocities  (hat  can  only  be 

zero-winds  simulation.  accounted  for  by  the  wind.  This  is  not  to  say  that 
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Figure  32.  Actual  and  cumulative  daily  predicted 
trajectories  for  buoy  1364  for  the  zero-winds  si¬ 
mulation. 


currents  in  this  region  do  not  also  undergo  rapid 
temporal  variations.  It  is  probably  safe  to  assume, 
however,  that  the  response  of  the  currents  to  the 
wind  will  be  less  than  that  of  the  ice. 


the  midpoint  of  the  grid  is  very  well  predicted. 
This  feature  was  somewhat  apparent  in  the  60-day 
thickness  field  produced  by  the  standard  run  (Fig. 
Sg),  but  was  partially  obscured  by  the  new  ice 
growth.  Flere  it  is  clear  that  (his  feature  in  the  ob¬ 
served  ice  edge  is  probably  a  result  of  ice 
dynamics. 

That  too  little  ice  is  present  in  the  north  is  al¬ 
most  certainly  a  result  of  no  ice  growth.  The  ex¬ 
cess  of  ice  in  the  south  is  an  unresolved  problem, 
but  two  possibilities  exist.  First,  the  simulated  ice 
dynamics  may  not  be  adequately  reproducing  the 
actual  ice  dynamics  in  (his  area,  due  to  improper 
winds,  currents,  or  possibly  ice  rheology.  The 
other  possibility  is  that  the  ocean  significantly  ab¬ 
lated  the  ice  in  this  region  during  this  time  period. 

The  thickness  buildups  along  the  coast  are  more 
numerous  than  those  of  the  standard  run,  presum¬ 
ably  due  to  slightly  higher  velocities.  This  figure 
also  shows  many  smaller  areas  of  lower  concen¬ 
tration  and  thickness,  however,  because  new  ice 
growth  was  not  allowed  to  proceed.  Offshore 
velocities  and  subsequent  mass  divergence  in  the 
lee  of  boundary  promontories  appear  to  be  re¬ 
sponsible  for  the  lower  concentrations  along  the 
coast.  Away  from  the  coast,  excessive  advection 
out  of  cells  creates  the  divergent  areas. 

In  light  of  the  superior  ice  extent  prediction,  it  is 
meaningful  to  again  provide  a  table  of  total  ice 
coverage  as  has  been  done  previously.  Table  3 
shows  the  total  area  coverage  for  the  dynamics 
simulation  at  10-day  intervals,  the  percentage 
change  during  intervals,  and  the  percent  differ¬ 
ence  between  this  simulation’s  predictions,  the  ob¬ 
servations  and  the  standard  run. 


Dynamics  simulation 

When  it  became  apparent  that  the  thermody¬ 
namics  was  dominating  the  ice  extent  in  the  previ¬ 
ous  simulations,  it  was  decided  to  simulate  the 
60-day  period  without  thermodynamics.  The  idea 
here  was  to  see  if  the  ice  dynamics  alone  could  en¬ 
large  the  ice  extent,  and  if  significant  features 
caused  by  ice  dynamics  in  previous  simulations 
were  being  masked  by  the  overwhelming  ice 
growth.  For  this  simulation,  zero  growth  rates 
were  assigned.  Also,  the  geostrophic  current  field 
of  earlier  simulations  was  used. 

The  60-day  thickness  and  compactness  fields 
which  are  shown  in  Figure  33  are  rather  surpris¬ 
ing.  The  predicted  ice  edge  (0.2  compactness  con¬ 
tour)  matches  the  observed  ice  edge  better  than  in 
any  previous  simulation.  In  addition,  the  sharp 
break  of  the  actual  ice  edge  towards  the  coast  near 


Table  3.  Areas  of  ice  cover  for  Ike  dynamics  simu¬ 
lation  (10-day  intervals). 


Predicted 
ttO"  m’f 

%diff 

from  observed 

5k  diff  from 
std  run 

Initial  area 

1.80 

Area  day  10 

1.65 

-25.3 

-12.2 

Change  (V») 

-8.3 

Area  day  20 

1.83 

-18.3 

-48.4 

Change  (%> 

10.9 

Area  day  30 

1.85 

-10.2 

-45.2 

Change  (4%) 

l.l 

Area  day  40 

1.98 

-15.0 

-47.8 

Change  (89) 

7.0 

Area  day  30 

2.26 

-8.9 

-49.7 

Change  <49) 

14.1 

Area  day  60 

2.37 

-24.2 

-59.4 

Change  (49) 

4.9 
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Figure  33.  Sixty -day  thickness  and  compactness  fields  for  the  dynamics  simulation.  Dashed  line  is  observed  ice 
edge  position  for  2  December  1979  (NPOC  1979). 


Generally,  differences  between  the  dynamics 
simulation  predicted  coverage  and  the  observed 
coverage  are  far  less  than  those  of  the  standard 
simulation.  As  expected  from  previous  figures,  the 
extent  of  ice  is  also  far  less  than  that  simulated  by 
the  standard  run.  The  less-than-observed  extents 
predicted  here  do  indicate,  however,  that  ice 
growth  is  necessary  for  a  reasonable  simulation. 
This  is  made  quite  obvious  by  the  particularly 
large  differences  between  the  simulated  and  ob¬ 
served  extents  on  day  60.  Likewise,  the  large  ex¬ 
pansion  of  the  ice  cover  that  actually  took  place 
between  days  50  and  60  (26%  increase)  appears  to 
be  primarily  due  to  ice  growth.  The  thermody¬ 
namic  simulation  predicted  such  an  increase 
(30%),  but  the  ice  cover  in  that  simulation  was  al¬ 
ready  excessive  on  day  50  and  the  increase  extend¬ 
ed  the  cover  too  far  eastward. 

The  implication  here  has  been  made  previously. 
This  simulation  verifies  that  both  dynamics  and 
thermodynamics  are  important  to  obtain  reasona¬ 
ble  ice  extents  in  this  region.  The  thermodynamics 
used  for  these  simulations  needs  considerable  im¬ 


provement,  however.  It  seems  that  the  major  im¬ 
provement  to  be  made  will  be  an  adequate  specifi¬ 
cation  of  oceanic  heat  flux.  That  the  thermody¬ 
namic  code  is  adequate  for  the  Arctic  Basin  has 
been  shown  by  Hibler  (1980b).  In  that  study,  oce¬ 
anic  boundary  layer  heat  storage  was  allowed  and 
the  simulation  provided  very  reasonable  results.  In 
the  East  Greenland  area,  perhaps  both  boundary 
layer  heat  storage  and  the  advection  of  warmer 
waters  from  the  south  need  to  be  parameterized 
for  an  adequate  thermodynamic  model. 

The  60-day  averaged  velocities,  shown  in  Figure 
34,  look  similar  to  those  of  the  standard  run  (Fig. 
7g)  except  for  areas  near  the  eastern  side  of  the 
high  velocity  stream.  In  the  region  of  the  sharp 
break  toward  the  coast,  an  onshore  velocity  com¬ 
ponent  exists.  In  previous  simulations,  higher 
thicknesses  in  this  region  presumably  prevented 
this  onshore  component.  Velocities  within  the 
stream  itself  also  appear  to  be  slightly  higher,  but 
these  differences  are  difficult  to  discern  from  these 
figures. 

The  predicted  u  and  v  velocity  component 
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Figure  34.  Sixty-day  averaged  velocity  field  for  the 
dynamics  simulation. 


Figure  35.  Actual  and  cumulative  daily  predicted 
trajectories  for  buoy  1564  for  the  dynamics  simu¬ 
lation. 


means  for  buoy  1564, 0.08  and  -0.18  m  s*1,  exceed 
those  in  any  previous  simulation  as  far  as  compar¬ 
ison  to  the  observed  component  means  (0.08  and 
-0.19  m  s'1)  is  concerned.  In  addition,  the  RMS  er¬ 
rors  for  the  velocity  components  are  quite  reason¬ 
able,  being  0.13  and  0.14  m  s'1.  The  u  component 
error  is  0.01  m  s'1  higher  than  that  for  the  stan¬ 
dard  run  while  the  v  component  error  is  the  same 
as  in  the  standard  simulation.  In  the  correlation  of 
predicted  velocities  with  those  of  buoy  1564,  the  u 
and  v  coefficients  are  0.49  and  0.53.  These  coeffi¬ 
cients  represent  a  slightly  higher  value  for  the  u 
component  and  a  lower  value  for  the  v  component 
than  in  the  standard  run  (0.48,  0.57)  but,  in 
general,  are  in  the  same  range  as  the  other  simula¬ 
tions.  The  velocity  differences  become  more  ap¬ 
parent  in  the  cumulative  daily  predicted  trajectory 
of  this  buoy,  shown  in  Figure  35.  This  trajectory  is 
better  than  that  obtained  in  any  other  simulation 
with  respect  to  final  position.  This  is  probably  ac¬ 
counted  for  by  slightly  higher  velocities  overall. 

The  dynamics  simulation  predicts  a  14%  larger 
volume  of  northern  inflow  ice  than  the  standard 


run,  5.39x10"  m1.  This  is  no  surprise  in  that 
larger  velocities  would  be  expected  in  the  inflow 
region  because  of  lower  strength  ice.  The  lower 
strengths  presumably  result  from  lower  thickness¬ 
es  due  to  the  lack  of  ice  growth.  The  southern 
outflow,  on  the  other  hand,  differs  considerably. 
The  dynamics  simulation  produces  0.32x10"  m’ 
while  the  standard  run  predicts  2'/2  times  this 
amount.  This  is  obviously  a  result  of  the  lack  of 
thermodynamics,  whose  effect  is  twofold.  Ini¬ 
tially,  the  ice  growth  acts  to  ‘Till”  the  region  with 
ice  quite  rapidly  (too  rapidly)  and  outflow  for  the 
standard  run  can  begin  at  an  earlier  time.  In  addi¬ 
tion,  the  continuing  growth  of  ice  makes  consider¬ 
ably  more  ice  available  for  outflow  during  the  en¬ 
tire  simulation  period. 


SUMMARY  AND 
CONCLUDING  REMARKS 

Pertinent  mass  balance  results  from  the  various 
simulations  are  presented  in  Table  4.  Velocity 
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Table  4.  Model  sensitivity  test  mass  balance  results. 


Northern 
inflow 
HO"  in') 

Southern 
outflow 
110"  m') 

Total 

urowth 

1 10"  in') 

Total  ice 

increase 

HO"  m  i 

60-day  areal 
coverage 

HO  m') 

Standard 

4.73 

0.83 

3.95 

7.85 

5.84 

Thermodynamics 

0.0 

0.0 

2.57 

2.57 

5.56 

Zero  strength 

6.18 

-1.1 1 

155.56 

162.85 

4.47 

Zero  inflow 

0.0 

0.83 

5.62 

4.77 

5.56 

Zero  currents 

3.85 

0.81 

4.11 

7.16 

5.78 

Modified  currents 

7.41 

2.94 

5.61 

10.08 

5.88 

Zero  winds 

2.02 

0.45 

2.74 

4.31 

6.15 

Dynamics 

5.39 

0.32 

0.0 

5.07 

2.37 

Table  5.  Simulated  velocity  comparisons  with  buoy  1564. 


Predicted  RMS  error  Correlation 

mean  ♦  (m  s')  (m  s')  coefficient 


u 

V 

14 

V 

14 

V 

Standard 

0.08 

-0.17 

0.12 

0.14 

0.48 

0.57 

Zero  strength 

-0.03 

-0.17 

0.16 

0.16 

0.45 

0.45 

Zero  inflow 

0.06 

-0.16 

0.14 

0.15 

0.48 

0.52 

Zero  currents 

0.07 

-0.15 

0.12 

0.14 

0.47 

0.54 

Modified  currents 

0.13 

-0.25 

0.20 

0.21 

0.42 

0.50 

Zero  winds 

0.06 

-0.07 

0.09 

0.17 

0.03 

0.29 

Dynamics 

0.08 

-0.18 

0.13 

0.14 

0.49 

0.53 

*  The  observed  u  and  v  means  are  0.08  and  -0.19  m  s'1. 


comparison  statistics  for  each  simulation  are  sum¬ 
marized  in  Table  5.  Most  of  these  have  been  previ¬ 
ously  discussed  but  are  repeated  here  for  the  sake 
of  ease  of  comparison. 

By  far  the  most  radical  statistics  are  yielded  by 
the  zero  strength  simulation.  Total  ice  growth,  as 
previously  mentioned,  is  two  orders  of  magnitude 
higher  than  that  predicted  by  any  of  the  other  sim¬ 
ulations.  This  is  presumably  produced  by  the  con¬ 
tinuing  cycle  of  rapid  advection  and  subsequent 
thin  ice  growth  in  specific  areas.  From  a  different 
perspective,  the  free  drift  condition  causes  the  ice 
to  respond  nearly  instantaneously  to  rapidly  vary¬ 
ing  winds,  regardless  of  the  ice  thickness.  In  other 
simulations,  this  effect  is  diminished  by  allowing 
the  ice  to  have  strength.  This  results  in  a  more 
uniform  velocity  field.  While  areas  of  divergence 
can  still  occur,  it  is  only  after  periods  of  sustained 
wind  forcing  that  ice  is  moved  into  areas  of  less 
strength.  Another  surprising  result  of  the  zero 
strength  simulation  is  that  net  inflow  took  place  in 
the  southern  boundary  region,  although  this  is  not 
apparent  from  the  averaged  velocity  field  (Fig. 
20).  A  detailed  examination  of  the  results  resolved 
this  problem.  During  the  latter  part  of  the  simula¬ 


tion,  reasonably  thick  ice  had  accumulated  adja¬ 
cent  to  the  coast  in  this  boundary  region.  In  addi¬ 
tion,  a  slight  northward  component  of  velocity 
was  evident  for  several  grid  points  in  that  area.  As 
a  result,  this  small  northward  drift  moved  the  very 
thick  ice  back  into  the  grid  and  a  net  inflow  in  this 
region  was  created.  The  essence  of  this  simulation 
is  that  the  zero  strength  condition  allows  the  forc¬ 
ing  fields  to  move  unreasonably  large  amounts  of 
ice.  It  is  apparent  that  an  ice  strength  that  is  at 
least  partially  dependent  on  thickness  is  necessary 
to  obtain  reasonable  ice  thicknesses  and  compact¬ 
nesses. 

The  thermodynamics  simulation  made  it  clear 
that  ice  growth  alone  was  responsible  for  creating 
too  large  an  area  of  ice  cover.  It  became  apparent 
in  this  simulation  that  the  oceanic  heat  flux  needs 
to  be  considered  to  properly  model  ice  growth  in 
this  region.  Comparison  of  the  total  growth  in  this 
simulation  with  that  produced  in  the  standard  run 
also  shows  that  ice  growth  is  further  stimulated  by 
ice  dynamics.  This  fact  may  have  a  significant  ef¬ 
fect  on  air-sea  energy  exchange  in  this  region.  For 
this  reason,  and  for  the  obvious  reason  that  ice 
transport  and  drift  cannot  be  predicted  with  a 


pure  thermodynamic  ice  model,  the  inclusion  of 
ice  dynamics  in  any  modeling  effort  for  this  region 
is  deemed  necessary. 

The  zero  inflow  simulation  emphasizes  the  im¬ 
portance  of  ice  transport  from  the  Arctic  Basin  in¬ 
to  the  East  Greenland  Sea.  Although  the  areal 
coverage  is  not  significantly  different  from  that 
predicted  by  the  standard  simulation,  the  thick¬ 
ness  and  compactness  fields  (Fig.  22)  and  the  total 
ice  volume  are  quite  different.  No  impact  upon  the 
volume  of  southern  outflow  is  noted,  but  a  longer 
simulation  would  likely  show  discernible  differ¬ 
ences.  Reasonable  ice  velocities  are  maintained  in 
this  simulation,  and  as  a  result,  total  growth  is 
larger  than  that  of  the  standard  run.  This  is  due 
primarily  to  the  rapid  advection  and  subsequent 
growth  in  the  northern  sector  where  velocities  are 
larger.  The  salient  point  of  this  simulation  is  that  a 
reasonable  thickness  regime  cannot  be  maintained 
without  ice  inflow  which,  when  corrected  for 
outflow,  supplied  approximately  half  the  total  ice 
volume  increase  during  this  period  (according  to 
the  standard  run). 

The  zero  currents  simulation  sheds  light  on  the 
fact  that  the  geostrophic  currents  used  in  other 
simulations  contribute  little  to  the  ice  dynamics 
except  to  the  volume  of  inflow  ice.  In  view  of  pre¬ 
vious  investigations  reviewed  earlier,  it  is  sus¬ 
pected  that  the  geostrophic  current  velocities  are 
too  small.  On  the  other  hand,  using  the  average  ice 
velocity  field  as  a  current  field  creates  excessive  ice 
velocities.  These  large  velocities  greatly  increase 
inflow,  outflow  and  ice  growth  to  levels  that  are 
probably  also  unreasonable.  What  these  two  simu¬ 
lations  tend  to  point  out  is  that  if  all  other  model 
parameters  are  reasonable,  then  the  currents  in 
this  region  are  neither  purely  geostrophic  (also  as¬ 
suming  the  geostrophic  currents  here  are  reason¬ 
able)  nor  totally  ice-driven.  This  is  not  unex¬ 
pected,  and  future  work  at  some  point  should  ad¬ 
dress  a  coupled  ice-ocean  model. 

The  zero  winds  simulation,  which  uses  the  mod¬ 
ified  current  field  to  create  an  extreme  case,  clear¬ 
ly  shows  that  the  winds  used  in  this  study  were  re¬ 
sponsible  for  most  of  the  ice  transport  during  this 
simulated  time  period.  Even  with  the  excessive 
currents,  the  total  transport  (inflow  and  outflow) 
is  about  half  that  simulated  by  winds  alone.  Daily 
winds  are  also  necessary  to  predict  reasonably  ac¬ 
curate  velocities,  as  noted  from  the  large  daily  va¬ 
riations  in  actual  buoy  velocities.  Although  it  ap¬ 
pears  that  ice  dynamics  can  be  reasonably  simu¬ 
lated  without  currents  for  short  time  periods, 
these  results  imply  that  they  cannot  be  properly 
simulated  without  winds.  As  mentioned  earlier. 


the  wind  velocities  may  be  excessive  near  the  coast 
and  too  small  in  the  ice  edge  region.  More  detailed 
analysis  of  the  wind  fields  will  be  the  subject  of  fu¬ 
ture  investigations. 

That  ice  growth  has  been  excessive  in  previous 
simulations  is  verified  by  the  dynamics-only  case. 
In  this  simulation  the  best  agreement  between  the 
observed  and  predicted  ice  extent  is  obtained. 
From  the  60-day  thickness  and  compactness  con¬ 
tour  plot  (Fig.  31)  it  is  obvious,  however,  that 
growth  is  necessary  in  the  north,  and  more  abla¬ 
tion  is  needed  in  the  southern  region  to  obtain  a 
better  predicted  areal  extent.  The  volume  of  in¬ 
flow  ice  in  this  simulation  appears  excessive.  As 
was  apparent  from  the  standard  simulation, 
growth  in  the  northern  sector  sufficiently  in¬ 
creased  thickness,  which  decreased  velocities  and 
suppressed  inflow.  The  southern  outflow  volume 
here  is  approximately  40%  of  that  in  the  standard 
run.  It  is  suspected  that  this  value  is  the  more  rea¬ 
sonable,  due  primarily  to  excessive  ice  extent  in 
the  south  when  growth  is  allowed  to  proceed. 

The  problems  with  the  model  results,  in  particu¬ 
lar  referring  to  the  standard  run,  have  yet  to  be  re¬ 
solved.  As  noted  on  several  occasions,  these  major 
problems  are  excessive  ice  growth  and  improper 
ice  velocities  in  the  vicinity  of  the  ice  edge.  Future 
work  will  focus  on  the  resolution  of  these  prob¬ 
lems  within  guidelines  previously  mentioned.  In 
addition,  simulations  will  be  carried  out  for  differ¬ 
ent  seasons  and  hopefully  for  longer  periods 
(90-120  days). 

In  spite  of  these  problems,  however,  this  study 
has  shed  light  on  several  key  issues  concerning 
modeling  studies  in  the  East  Greenland  area. 
Without  hesitation,  the  most  important  is  that  a 
sea  ice  model  which  utilizes  a  viscous-plastic  con¬ 
stitutive  law  as  developed  by  Hibler  (1979)  seems 
to  provide  reasonable  results  over  most  of  the  re¬ 
gion.  Furthermore,  this  investigation  points  out 
the  necessity  of  using  a  coupled  dynamic-thermo¬ 
dynamic  model  to  properly  model  this  region.  The 
importance  of  including  ice  dynamics  in  studies  of 
air-sea  energy  exchange  has  previously  been  em¬ 
phasized.  Allowing  the  ice  to  have  strength  and  to 
interact  with  itself  is  a  necessity  that  has  also  been 
clarified  by  this  study.  In  addition,  it  has  been 
shown  that  winds,  currents  and  ice  import  from 
the  Arctic  Basin  all  contributed  significantly  to  the 
ice  balance,  even  during  (his  short  study  period. 

Detailed  model  refinement  will  be  difficult  with¬ 
out  a  significant  increase  in  the  amount  and  qual¬ 
ity  of  the  observational  data.  In  particular,  future 
studies  will  require  more  ice  drift  data  and  more 
detailed  thickness  and  concentration  information. 


The  relatively  sparse  data  set  used  in  this  study  has 
been  sufficient  to  point  out  certain  model  deficien¬ 
cies  and  to  draw  general  conclusions.  The  model 
can  also  be  further  refined  to  a  limited  degree  (i.c. 
thermodynamics  and  possible  ice  velocities)  with 
the  currently  available  data.  However,  the 
refinements  necessary  to  “tune”  the  model  for 
operational  forecasting  use  will  certainly  require 
more  detailed  data  during  all  seasons. 
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